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Abstract 


This  document  describes  in  detail  a  software  tool  for  manipulating  data  files.  The 
Surveillance  Acoustics  section  at  Defence  Research  Establishment  Atlantic  has  acquired 
VAX  computers  over  the  last  few  years,  and  analysis  tasks  which  were  formerly  done  on 
PDP-1 1  computers  are  now  being  moved  to  the  VAXen.  PDP-1  Is  are  still  used  in  the  at- 
sea  data  collection  role,  so  some  means  is  necessary  of  transferring  the  data  files  thus 
produced  to  the  VAXen  for  signal  processing  and  analysis.  PDP-1 1  data  files  are  typically 
located  on  9-track  magnetic  tape,  so  one  method  of  transferring  the  data  would  be  to  read 
PDP-11  tapes  on  the  VAXen.  The  software  tool  described  here  (a  program  named 
TRANSFER)  was  written,  in  part,  to  perform  this  data  transfer  chore,  taking  into  account 
the  special  formats  and  header  information  in  the  files  produced  by  the  PDP-1  Is. 
Manipulation  of  data  files  already  residing  on  a  VAX  is  also  possible  using  TRANSFER. 
The  program  is  versatile,  allowing  the  user  to  choose  channels  and  data  segments  to  be 
transferred  between  files  with  a  high  degree  of  freedom. 


Sonunaire 


Le  present  article  deem  en  detail  un  outil  logiciel  permettant  la  manipulation  de 
fichiers  de  donnees.  La  section  de  l'acoustique  de  surveillance  du  Centre  de  recherches 
pour  la  defense.  Atlantique.  sest  dote  d  ordinateurs  VAX  au  cours  des  demieres  annees  et 
les  travaux  danalyse  effectues  jusqua  present  sur  des  PDP-11  le  sont  maintenant  au 
rnoyen  du  VAXen.  Toutefois,  les  PDP-11  sont  encore  utilises  pour  la  cueillette  des 
donnees  en  mer:  les  travaux  d  analyse  et  de  traitment  des  donnees  necessitent  done  une 
methode  de  transfert  sur  le  VAXen  des  fichiers  produits  pendant  la  cueillette  des  donnees. 
Puisque  les  bandes  magnetiques  de  9  pistes  constituent  le  support  de  me  moire  typique  des 
fichiers  de  donnees  PDP- 1 1 .  on  peut  envisager  la  procedure  suivante  comme  methode  de 
transfert  des  donnees:  lecture  des  donnees  memorisees  sur  les  bandes  PDP-11,  puis 
transfert  au  VAXen.  L’outil  logiciel  deent  dans  le  present  texte  tun  proeramme  appele 
TRANSFER!  a  ete  ecrit  en  partie  en  vue  d  assurer  cette  tache  de  transfert  des  donnees.  en 
fonction  des  donnees  d'entete  et  des  formats  speeiaux  des  fichiers  produits  par  le  PDP- 1 1 
TRANSFER  permet  aussi  la  manipulation  des  fichiers  de  donnees  d£ ja  memorises  dans  un 
ordinateur  \' AX.  Ce  logiciel  est  ires  polyvalent  et  permet  a  1  utilisateur  de  clioisir  avec  un 
crand  degre  de  liberie  les  pistes  et  les  segments  de  donnees  a  transferer  d  un  fichier  a 
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I.  INTRODUCTION 


Until  recently,  data  collected  at  sea  by  DREA  scientists  were  analysed  almost 
exclusively  on  PDP-1 1  /34s.  When  the  Surveillance  Acoustics  section  began  to  purchase 
VAX  computers,  the  situation  changed  drastically.  Most  analysis  tools  which  existed  on 
the  1  l/34s  are  now  available  on  the  VAXen  (the  accepted  plural  form  for  VAX)  in  a  similar 
or  more  powerful  form.  New  tools  (such  as  a  suite  of  shot-analysis  programs)  have  also 
been  appearing.  The  multi-user  nature  of  the  VAXen  has  made  these  software  tools 
available  to  a  wider  range  of  users,  and  consequently  more  analysis  is  being  performed. 

PDP-1  Is  have  been  used  in  the  at-sea  data  collection  role  at  DREA  for 
approximately  10  years,  and  perform  their  job  very  efficiently.  Real  time  data  collection 
programs  write  data  (generally  to  9  track  magnetic  tape)  in  a  format  which  has  been  well 
tuned  to  the  needs  of  DREA  scientists  over  the  lifetime  of  the  PDP-1  Is.  Naturally  enough, 
this  format  is  somewhat  foreign  to  the  VAX  computer,  and  data  tapes  produced  on  the 
PDP- 1 1  s  cannot  be  read  directly  on  the  VAX  without  a  software  interface.  To  allow  more 
users  to  access  and  analyze  raw  data,  such  a  software  interface  has  been  written,  and  this 
note  descri  bes  that  tool,  a  program  called  TRANSFER. 

TRANSFER  has  been  written  to  be  as  general  as  possible,  leaving  many  options 
open  to  the  user.  Most  files  written  to  tape  by  the  PDP-11  data  acquisition/analysis 
programs  can  be  moved  from  mag-tape  to  the  VAX  using  this  program.  The  formats 
accepted  are  ’  .DAT"  (time  senes  data  files),  ".FTR”  (Fourier  coefficient  files ;,  and 
•  PWR"  (power  spectrum  files).  Other  features  such  as  transferring  a  segment  of  an  input 
file  defined  by  start  and  stop  times,  or  transferring  only  a  subset  of  the  total  number  of 
channels  in  the  input  file  are  also  available.  Disk  space  is  at  a  premium  on  the  VAXen  (as  it 
is  on  all  computers),  so  the  latter  feature  is  an  important  one.  It  allows  disk  space  to  be 
conserved  if  a  user  w  ishes  to  analyze  only  a  few  of  the  channels  of  input  data  available. 

The  default  format  for  VAX  disk  files  created  by  TRANSFER  is  binary  direct 
access.  This  format  allows  random  access  to  any  block  in  the  file,  and  is  a  relatively 
compact  storage  format.  Most  of  the  analysis  programs  now  present  on  the  VAXen  accept 
input  files  in  this  format.  An  older  format  (READRT)  is  still  used  by  some  analysis 
programs,  and  TRANSFER  will  use  this  as  its  output  format  if  requested. 

Size  reduction  for  disk  files  already  on  the  VAX  is  desirable  in  many 
circumstances,  so  TRANSFER  will  also  perform  VAX  disk  file  to  VAX  disk  file  transfers. 
The  same  options  are  offered  in  the  disk-to-disk  mode  as  in  the  tape-io-disk  mode. 

The  next  section  of  this  note  describes  the  formats  of  the  tape  files  TRANSFER  will 
accept  from  PDF’-]  Is.  and  gives  more  details  on  the  VAX  and  READRT  disk  file  formats. 
Follow  ing  that,  a  detailed  description  of  program  implementation  will  be  given,  including  a 
discussion  of  all  the  options  available.  An  example  of  program  use  is  then  presented  to 
give  the  reader  some  feeling  of  how  a  terminal  session  proceeds  and  finallv,  some 
possibilities  for  future  developments  of  the  program  are  proposed. 


2.  FILE  FORM  ATS 

TRANSFER  was  originally  written  to  accept  P-track  magnetic  tapes  written  bv  the 
DREA  PDP-1  1  data  recording  programs,  and  transfer  them  to  VAX  files.  To  make  tho 
divument  self-contained,  a  brief  description  of  the  format  of  the  various  types  of  tape  input 


files  will  be  presented  here.  A  fuller  description  of  the  file  formats  and  the  philosophy 
behind  their  structure  is  contained  in  a  Technical  Communication  by  D.  Caldwell  [1].  Disk 
file  types  which  are  compatible  with  TRANSFER  will  also  be  described. 

2.1  Time  Series  (.DAT)  File  Format 

Time  series  data  (.DAT)  files  are  the  most  common  type  of  input  file  used  with  the 
TRANSFER  program.  Time  series  data  are  recorded  on  magnetic  tape  on  the  PDP-1  Is  in 
this  format,  and  since  it  is  the  intention  to  use  the  VAX  for  most  analysis,  raw  input  data 
will  be  moved  to  the  VAX  via  this  file  type. 

Each  tape  file  begins  with  a  512  byte  header  which  describes  the  physical 
parameters  of  the  file  such  as  record  length,  sampling  frequency,  etc.  The  header  is 
divided  into  four  blocks  -  the  first  32  bytes  form  an  integer  block  (2  bytes  per  integer),  the 
next  32  bytes  form  a  floating  point  block  (4  bytes  per  floating-point  number),  the  next  128 
bytes  form  a  byte  block  (1  byte  per  entry)  and  the  remaining  320  bytes  form  an  ASCII 
block  (1  character  per  byte).  Table  I  illustrates  these  blocks  and  gives  a  brief  description 
of  the  contents  of  each  location  in  the  header  block.  A  more  detailed  description  of  the 
meaning  of  the  header  block  contents  can  be  found  in  [1]. 

The  data  portion  of  a  typical  time  series  file  is  diagrammed  in  Figure  1.  In  the 
sample  file  shown,  there  are  m  time  samples  for  each  of  the  n  input  channels.  The  values 
are  multiplexed  so  that  the  first  time  sample  for  each  channel  appears  in  sequence,  followed 
by  the  second  time  sample  for  each  channel,  etc.  The  DREA  header  actually  allows  the 
time  series  data  to  be  written  in  other  formats,  but  the  one  shown  here  is  used  almost 
exclusively. 

22  Fourier  Coefficient  (.FTR)  File  Format 

This  type  of  file  is  used  by  analysis  programs  which  require  Fourier  coefficients  but 
do  not  contain  an  FFT  module  of  their  own.  Programs  for  performing  interference 
cancellation  which  require  ".FTR"  files  currently  exist  on  the  VAX.  A  program  (called 
SAFTR)  [2]  can  be  used  to  produce  ".FTR"  files  from  ".DAT"  files  on  the  VAX. 
TRANSFER  will  work  with  the  ".FTR"  format,  but  the  ".DAT"  format  is  more  likely  to  be 
encountered. 

As  in  the  ".DAT"  format,  each  ".FTR"  file  begins  with  a  512  byte  header  [1],  Table 
II  gives  a  brief  description  of  the  contents  of  each  location  in  the  ".FTR”  header  block. 
Figure  2  shows  the  contents  of  a  typical  ".FTR"  file.  In  the  sample  file  shown,  the  data  are 
multiplexed  in  a  different  manner  from  that  in  a  ".DAT"  file.  Here,  all  Fourier  coefficients 
from  each  transform  of  each  channel  are  kept  together;  that  is,  blocks  of  data  for  each 
channel  rather  than  single  samples  are  multiplexed. 

23  Power  Spectrum  (.PWR)  File  Format 

A  power  spectrum  tape  file  is  likely  to  be  transferred  from  tape  only  if  time  series 
analysis  was  performed  on  a  PDP-1 1.  SEQFFT  is  the  most  widely  used  spectral  analysis 
program  on  the  PDP- 1  Is,  and  the  output  of  that  program  conforms  to  the  ".PWR"  format. 
As  analysis  effort  moves  to  the  VAX  from  the  11s,  transfer  of  this  file  type  is  likely  to 
become  less  common. 

The  ".PWR"  file  begins  with  a  512  byte  header  [1].  TABLE  III  gives  a  brief 
description  of  the  contents  of  each  location  in  the  ".PWR"  header  block.  Figure  3  shows 


INTEGER  BLOCK 


I LABEL ( 1 )  -  Block  size 
I LABEL (2)  —  Record  Size 
ILABEL(3)  -  #  of  records 
ILABEL(4)  -  Repetition  rate 
ILABEL (5)  -  Number  type 
ILABEL ( €)  —  Bytes  per  number 
ILABEL(7)  —  #  of  channels 
ILABEL (8)  —  Multiplex  length 
ILABEL (9)  -  #  accumulations 
ILABEL (10)-  X-axis 
ILABEL (11)-  Y-axis 
ILABEL (12)-  Sequence  # 

ILABEL (13)-  Block  Scaling 
ILABEL (14)-  Spare 
ILABEL (15)-  Spare 
ILABEL (16)-  Spare 


#  of  16  bit  words/physical  block 

#  of  words  per  logical  record 

#  of  data  records  in  file 

#  of  records  per  repetition  cycle 
Int-1, Flt-2 , Cmplxl-lloct, CF-12oct 
eg.:  1-2 , F-4 , CI-4 , CF-8 

Must  divide  evenly  into  ILABEL (2) 

Word-1,  Record-Record  size 

Normally  1 

Time-1,  Frequency-2 

Linear-1, Square-2,  Log— 4 

User  assigned,  usually  increments 

Power  of  2  scaling  factor 


FLOATING-POINT  BLOCK 


FLABEL (1) 
FLABEL (2) 
FLABEL (3) 
FLABEL (4) 
FLABEL (5) 
FLABEL ( 6 ) 
FLABEL ( 7 ) 
FLABEL ( 8 ) 


Sampling  freq(Hz)  -ve  means  heterodyned 
Heterodyning  freq.  -ve  means  real  heterodyned 
Reference  level  Calibration  factor 
Max.  magnitude 

Gain  correction  1.00343332  for  power  of  2 
Spare 
Spare 
Spare 


BLABEL (1) 
BLABEL (2 ) 
BLABEL (3) 
BLABEL (4) 


Channel  # 
Gain  <dB) 
Channel  # 
Gain  (dB) 


BYTE  BLOCK 

for  channel  in  BLABEL (1) 
for  channel  in  BLABEL (3) 


BLABEL (127)-  Channel  # 

BLABEL (128)-  Gain  (dB)  for  channel  in  BLABEL (127) 

ASCH  BLOCK 

ALABEL ( 1 )  -  First  character  of  ASCII  label  block 


ALABEL ( 320) —  Last  character  of  ASCII  label  block 


Table  I:  Header  contents  of  a  Time  Series  (.DAT)  file 


sample  1 

sample  2 

sample  3 

•  •  • 

sample  m 

Y  chan  1 

Y  chan  2 

Y  chan  3 

•  •  • 

Y  chan  n 

Figure  1:  Data  format  in  a  typical  Time  Series  (.DAT)  file 


INTEGER  BLOCK 

ILABEL(l)  -  Block  size  #  of  16  bit  words/physical  block 

ILABEL(2)  -  Record  Size  4  of  words  per  logical  record 

ILABEL(3)  -  #  of  records  4  of  data  records  in  file 

ILABEL(4)  -  Number  of  sequential  transforms 

ILABEL(5)  -  Number  type  Int-1 , Flt-2 , Cmplxl-lloct , CF-12oct 

ILABEL(6)  -  Bytes  per  number  eg.:  1-2 , F-4 , CI-4 , CF-8 
ILABEL(7)  -  #  of  channels  Must  divide  evenly  into  ILABEL(2) 

ILAB£L(8)  —  4  of  frequency  bins 
ILABEL(9)  -  4  accumulations  Normally  1 

ILABEL(IO)-  X-axis  Time-1,  Frequency-2 

ILABEL(ll)-  Y-axis  Linear-1, Square-2,  Log— 4 

ILABEL(12)—  Sequence  4  User  assigned,  usually  increments 

ILABEL(13)«  Block  scaling  Power  of  2  scaling  factor 

1LABEL(14)-  Window  type  None-1 , Hannina— 2 , Hamming-3 , Kaiser-4 

I LABEL ( 15 ) -  #  of  zeros 
ILABEL(16)-  #  of  points  of  overlap 


FLOATING-POINT  BLOCK 

FLABEL(l)  —  Start  frequency  of  first  bin  (Hz) 

FLABEL(2)  -  Heterodyning  freq.  -ve  means  real  heterodyned 
FLABEL(3)  -  Frequency  resolution  (Hz) 

FLABEL(4)  -  Max.  magnitude 
FLABEL (5)  -  Spare 

FLABEL ( 6 )  -  %  overlap 

FLABEL (7)  —  Time  interval  of  a  simgle  FFT  (hours) 

FLABEL (8)  -  Spare 


BYTE  BLOCK 


BLABEL ( 1 ) 
BLABEL ( 2 ) 
BLABEL (3) 
BLABEL ( 4 ) 


Channel  4 
Gain  (dB) 
Channel  4 
Gain  (dB) 


BLABEL (127)-  Channel  4 
BLABEL (128)-  Gain  (dB) 


ASCII  BLOCK 


ALABEL(l)  -  First  character  of  ASCII  label  block 


ALABEL(320)-  Last  character  of  ASCII  label  block 


Table  II:  Header  contents  of  a  Founer  Coefficient  (,FTR)  file 


DFTIchanl 

DFT1chan2 

DFTlchan3 

^1 

e-Q 

real  1 

imag  1 

real  2 

imag  2 

DFT2chan1  DFT2chan2 


real  p  imag  p 


Figure  2:  Data  format  in  a  typical  Fourier  Coefficient  (.FTR)  file 


INTEGER  BLOCK 

ILABEL(l) 

- 

Block  size 

•  of  16  bit  words /phys lea  1  block 

I LABEL (2 ) 

- 

Record  Size 

*  of  words  per  logical  record 

I LABEL ( 3 ) 

- 

#  of  records 

#  of  data  records  in  file 

I LABEL ( 4 ) 

- 

Number  of  sequential  spectral  estimates 

I LABE  L ( 5 ) 

- 

Number  type 

Int-1, Fit— 2 , Cmpixl-lloct,  CF-12cct 

I LABEL (6) 

- 

Bytes  per  number 

eg.:  1-2 , F-4 , Cl-4 , CF-6 

I LABEL ( 7 ) 

- 

*  of  channels 

I LABEL (6) 

♦  of  frequency  bins 

ILABEL  <  9 ) 

— 

*  accumulations 

Normally  1 

I LABEL (10) 

- 

X-axis 

Time-1,  Frequency-? 

ILABEL (11) 

- 

Y-axis 

Linear-i , Square— 2 , Loa— 4 

ILABEL ( 12 ) 

- 

Sequence  * 

User  assigned,  usually  increments 

ILABEL (13) 

- 

Spa  re 

ILABEL (14) 

- 

Window  type 

None— 1 , Hanning-2 , Hammina-3  ,  Ka 1 se  r-4 

ILABEL (15) 

— 

#  of  zeros 

ILABEL (16) 

- 

*  of  points  of  overlap 

FLOATING-POINT  BLOCK 

FLABEL (1) 

- 

Center  frequency 

of  first  bin  (Hz) 

FLABEL (2) 

- 

Heterodyning  freq 

-ve  means  real  heterodyned 

FLABEL (3) 

Frequency  resolution  (Hz) 

FLABEL (4) 

- 

Max.  magnitude 

FLABEL (5) 

- 

Spare 

FLABEL ( 6 ) 

- 

Time  interval  between  sequential  frames  (hrs . ) 

FLABEL (7) 

- 

%  overlap 

FLABEL (8) 

Spa  re 

BYTE  BLOCK 

BLABEL ( 1 ) 

- 

Channel  * 

BT  7. BEL  (2) 

- 

Normally  C . 

gain  already  compensated  for 

BLABEL '3) 

— 

Channel  # 

BLABEL ( * ) 

• 

Normally  0 

• 

BLABEL (127) 

■  Channel  * 

BLABEL (128) 

-  Normally  0 

ASCII  BLOCK 


ALABEL(l)  -  First  character  of  ASCII  label  block 


ALABEL(320)-  Last  character  of  ASCII  label  block 


Table  HI:  Header  contents  of  a  Power  Spectrum  (.PWR)  file 


Figure  3:  Data  format  in  a  typical  Power  Spectrum  (.PWR)  file 


the  contents  cf  a  typical  ".PWR"  file.  Multiplexing  of  the  data  is  similar  to  that  found  in 
.FTR"  files.  Here,  there  are  q  power  points  grouped  together  representing  the  power  in 
each  of  q  frequency  bins  for  each  channel.  Power  spectrum  levels  are  typically  stored  as 
decibels  (dB).  Once  again,  different  formats  for  the  data  storage  are  acceptable  as  long  as 
they  conform  to  what  is  described  in  the  file's  header,  but  the  format  of  Figure  3  is  by  far 
the  most  common. 

2.4  VAX  and  READRT  File  Formats 

Files  from  TRANSFER  are  stored  on  VAX  disks  in  one  of  two  formats.  The  most 
common  of  these  is  a  binary  direct  access  file.  This  format  is  compact  and  allow's  direct 
access  to  any  block  in  the  file. 

Data  files  are  typically  very  large,  often  taking  an  entire  2400  foot  reel  of  9  track 
magnetic  tape,  so  transferring  them  to  disk  causes  storage  problems  to  appear  quickly  if 
several  people  are  doing  analysis.  For  this  reason  the  disk  files  should  be  stored  in  as 
compact  a  form  as  possible. 

Random  access  to  any  pan  of  the  data  in  a  file  is  also  important.  A  scientist  doing 
time  series  analysis  may  be  interested  in  only  certain  segments  of  data  in  a  large  file,  so 
having  to  sequentially  access  each  record  to  get  to  the  desired  one  would  be  inefficient. 
After  analysis,  display  programs  (such  as  PLTPWR)  also  need  random  access  to  the  data 
so  that  any  segment  may  be  displayed  rapidly  and  in  any  order.  The  format  of  the  VAX 
files  permits  this. 

The  other  format  is  called  READRT  after  the  file  transfer  program  which  originally 
used  it.  This  format  is  used  by  some  of  the  older  analysis  programs,  and  results  in  "a 
sequential  unformatted  file.  It  does  not  take  a  great  deal  more  disk  space  than  the  VAX 
format,  but  the  advantages  of  random  access  to  data  are  lost.  It  is  not  recommended  that 
this  file  format  be  used  in  future  analysis  programs. 


3.  PRfXiRAM  implementation 

This  section  begins  with  a  description  of  the  program  structure  and  the  features 
available  in  TRANSFER.  Implementation  details  are  then  presented,  and  non-standard 
practises  are  discussed  more  fully. 

TRANSFER  was  written  in  FORTRAN  77,  but  is  not  easily  transportable  to 
computers  other  than  VAX  and  micro-VAX  models  made  by  Digital  Equipment 
Corporation  due  to  its  extensive  use  of  system  calls.  (System  calls  use  internals  of  the 
VMS  operating  system  directly  and  are  not  part  of  the  FORTRAN  77  standard.)  The 
system  calls  were  used  to  speed  up  tape  and  disk  access,  and  to  make  use  of  some  of  the 
powerful  capabilities  of  the  VAX-VMS  operating  system.  On  the  bright  side,  however,  the 
program  can  be  run  on  any  of  the  VAX  machines  from  DEC  without  modification. 
TRANSFER  can  be  used  from  any  ANSII  standard  computer  terminal,  but  works  best  w'ith 
a  VT1 00/200  series  terminal  or  emulator. 

Subroutines  from  many  sources  were  used  in  TRANSFER.  Asynchronous  disk 
input/output  routines  from  NRL  (the  Naval  Research  Laboratory)  in  Washington  D.C. 
proved  to  be  very  useful  in  this  implementation  [3],  Other  useful  routines  from  various 
groups  at  DREA  have  been  incorporated  to  avoid  duplication  of  programming  effort. 
Subroutines  obtained  from  outside  sources  will  be  noted  as  such  in  the  following 
discussion. 


3.1  Program  Structure 


TRANSFER  was  written  in  a  modular  format  in  order  to  facilitate  modifications  and 
additions.  'User-friendliness'  and  simplicity  of  use  were  major  considerations  in  program 
design.  Flexibility  is  a  keyword  for  TRANSFER  since  many  input  file  formats  must  be 
accessible  to  VAX  users,  and  the  data  in  those  files  should  be  easy  for  the  user  to 
manipulate.  Simplicity  of  use  and  a  high  degree  of  flexibility  are  not  always  compatible, 
but  the  attempt  has  been  made  to  achieve  both  objectives  with  TRANSFER. 

Even  the  most  efficiently  written  program  can  be  practically  useless  if  it  has  a  poor 
user  interface.  For  this  reason,  considerable  effort  was  put  into  making  the  user  interface 
of  TRANSFER  easy  to  understand  and  use.  The  terminal  input  session  has  been  separated 
into  related  modules  (for  example,  one  module  deals  with  defining  the  section  of  an  input 
file  to  be  transferred  into  the  output  file).  Each  module  is  presented  on  a  separate  screen  on 
the  user's  terminal,  and  a  heading  appears  at  the  top  of  the  screen  describing  the  purpose  of 
the  module.  Examples  of  this  will  be  given  in  a  later  section. 

Program  structure  is  outlined  in  the  flowchart  of  Figure  4.  There  are  three  basic 
segments  -  input,  processing,  and  output.  Within  these  segments,  subroutine  structures 
were  used  when  possible.  These  are  not  noted  in  the  chart,  but  more  detail  on  some  of 
them  will  be  given  in  a  later  section.  The  loop  structure  of  the  program  is  fairly  simple  at 
the  flowchart  level  but  became  rather  difficult  to  implement  because  of  the  differences  in 
input  file  structures  which  had  to  be  accommodated.  Most  of  the  options  available  to  the 
user  are  noted  in  the  chart,  and  will  be  discussed  in  detail  in  the  following  sub-section. 

3.2  Program  Features 

Input  files  for  TRANSFER  can  be  located  either  on  magnetic  tape  or  on  disk.  In 
the  case  of  tape,  ".DAT",  ".FTR"  and  ".PWR"  files  are  accepted,  while  only  VAX  format 
disk  input  files  are  accepted.  One  of  the  inconveniences  of  using  magnetic  tapes  is  the 
requirement  that  the  user  must  remember  to  allocate  the  tape  unit  and  mount  the  tape  (both 
VMS  commands).  TRANSFER  avoids  this  by  the  use  of  VMS  system  calls.  When  the 
user  specifies  a  tape  unit  (for  example  MSBO:)  as  part  of  an  input  filename,  TRANSFER 
tries  to  allocate  that  unit  and  then  mount  any  tape  found  on  it.  If  the  unit  is  already  assigned 
or  the  tape  cannot  be  mounted,  the  user  is  informed  and  program  execution  halts;  otherwise 
the  operation  is  transparent. 

The  program  will  not  stop  executing  if  it  encounters  a  parity  error  while  reading 
from  a  tape  file.  Parity  errors  on  ".DAT"  tape  files  are  a  distinct  possibility  due  to  the 
manner  in  which  the  PDP-1 1  data  collection  programs  operate.  High  speed  is  the  priority 
for  the  data  collection  programs,  so  no  error  checking  is  performed  while  writing  to  tape 
(resulting  in  the  possibility  of  parity  errors).  (Error  checking  uses  valuable  time  and  could 
cause  data  to  be  lost  in  some  cases.)  When  TRANSFER  hits  a  parity  error,  it  rewinds  the 
tape  to  the  last  good  record  of  data  and  substitutes  that  for  the  corrupted  data.  This  action 
could  bias  the  statistics  of  the  output  if  many  parity  errors  are  encountered,  but  was  deemed 
to  be  more  appropriate  than  destroying  the  time  synchronization  of  the  file  by  throwing 
away  data.  The  user  is  notified  each  time  a  tape  parity  error  is  encountered,  and  a  running 
count  of  the  errors  is  presented. 

TRANSFER  allows  wildcards  to  be  used  to  specify  a  family  of  input  files  with 
similar  names  (or  parts  of  names)  for  input  from  either  disk  or  tape.  A  short  description  of 
wildcarding  would  be  appropriate  for  the  uninitiated  -  it  hinges  upon  the  use  of  a  wildcard 
character  (in  this  case  '*').  Filenames  are  made  up  of  a  name  and  an  extension  separated 
by  a  period  (for  example  TEST.DAT).  The  family  of  files  with  the  name  TEST  and  any 
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Figure  4:  General  flowchart  for  TRANSFER 


extension  can  be  specified  by  typing  TEST.*.  Similarly  all  files  with  the  extension  ".DAT" 
can  be  specified  by  typing  *.DAT.  Typing  TES*.DAT  would  specify  all  files  having  a 
name  beginning  with  TES  and  having  die  extension  ".DAT".  When  TRANSFER 
encounters  a  wildcard  in  a  filename,  it  finds  all  files  which  satisfy  the  input  name  set  and 
processes  them  in  sequence.  A  user  can  specify  a  new  set  of  transfer  parameters  for  each 
input  file,  or  alternatively  can  set  up  parameters  for  only  the  first  file  and  use  those  same 
parameters  for  all  of  the  other  files. 

Tape  files  present  other  possibilities  as  well.  A  user  may  wish  to  transfer  a  number 
of  files  from  tape  without  having  to  specify  their  names  (for  example  the  second  through 
the  fifth  files  on  the  tape).  This  option  is  available  to  TRANSFER  users  through  an  input 
file  switch  (switches  are  only  available  when  using  input  files  on  magnetic  tape).  A  switch 
is  used  in  the  following  manner:  the  user  types  the  name  of  the  tape  unit  upon  which  the 
reel  is  mounted,  followed  by  a  slash  (/)  and  a  switch  parameter.  For  example,  to  transfer 
the  first  through  the  fifth  files  from  a  tape  on  unit  MSBO:,  the  user  would  type 
MSBO:/START=1/STOP=5  when  prompted  to  enter  an  input  filename.  Here  the  START 
switch  defines  the  file  on  tape  with  which  to  begin  the  transfer,  and  the  STOP  switch 
defines  the  number  of  the  final  file  which  is  to  be  transferred. 

There  is  only  one  other  switch  available;  that  is  the  /V  (or  verify)  switch.  When 
this  switch  is  included  after  an  input  tape  filename  (for  example  MSBO:TEST.DAT/V), 
TRANSFER  will  ask  the  user  for  verification  before  skipping  anv  file  which  it  encounters 
on  the  tape  (If  the  first  file  found  is  the  one  the  user  specified,  it  is  processed  without 
question).  This  switch  is  useful  if  the  user  wishes  to  process  the  first  file  on  a  tape  but 
doesn't  know  its  name.  In  that  case,  when  TRANSFER  asks  whether  it  should  skip  the 
file,  the  user  need  only  give  a  negative  reply  and  it  will  be  processed.  When  /V  is  not  used, 
the  entire  tape  will  be  scanned  for  a  filename  match  and  no  option  for  processing  non¬ 
matching  files  will  be  presented  to  the  user. 

The  default  filename  for  files  produced  by  TRANSFER  is  the  same  as  the  input 
filename  with  the  extension  .TFR  to  indicate  that  it  is  output  from  JRANSEER.  The  user 
can  supply  a  different  filename  and/or  extension  if  desired.  Output  file  format  (VAX  or 
READRT)  is  also  selectable. 

The  case  will  often  arise  where  only  a  small  portion  of  an  input  data  file  will  be 
used  for  analysis  or  display.  For  this  reason,  TRANSFER  allows  the  user  to  choose  a 
segment  of  the  input  file  for  transfer.  A  user-selected  segment  is  always  defined  by  its  start 
time.  The  time  at  which  recording  for  a  file  began  is  included  in  the  DREA  header.  This 
time  is  displayed  for  the  user,  who  can  then  specify  the  time  at  which  the  desired  segment 
of  data  begins.  TRANSFER  will  then  skip  into  the  file  to  the  desired  time  (a  rather 
complex  process,  since  the  time  and  number  of  records  to  be  skipped  depend  upon  the  type 
of  file  being  transferred1  ).  When  the  file  has  been  positioned  to  the  desired  start-time, 
the  user  is  given  the  option  of  specifying  the  segment  length  as  a  number  of  blocks  (512 
bvtes/block).  a  time  duration,  or  the  remainder  of  the  file.  Thus  a  user  has  great  control 
over  the  data  transferred  to  the  output  disk  file. 


Tor  example.  .FTR  and  .PWR  files  have  a  time  resolution  which  is  determined  by  the 
FFT  length  used  in  producing  the  file.  For  2kHz  samples  and  an  8K  FFT,  each  set  of 
Fourier  coefficients  covers  a  4  s  time  interval.  Due  to  the  construction  of  the  files,  a  time 
resolution  of  less  than  4  s  would  be  impossible  in  this  case.  In  this  situation,  the  actual 

start  time  is  the  accessible  time  closest  to  (but  not  less  than)  that  requested  by  the  user. 
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Full  control  over  the  channels  to  be  put  into  the  output  file  is  also  a  necessity. 

*  Channel  numbers  which  were  used  in  data  recording  are  entered  in  the  byte  section  of  the 

DREA  header  and  so  are  available  during  TRANSFER  operation.  The  user  can  select  any 
subset  (or  all)  of  the  channels  in  the  original  input  file  for  transfer  into  a  VAX  file.  The 
header  of  the  new  file  thus  created  will  be  modified  to  include  only  channel  numbers  of 
those  channels  presently  in  the  file.  This  capability  is  useful  in  reducing  the  amount  of  data 
stored  on  disk  since  non-acoustic  channels  or  channels  known  to  contain  corrupted  data 
need  not  be  transferred  to  disk. 

During  operation,  TRANSFER  provides  feedback  to  a  status  screen  on  the  user's 
terminal.  This  screen  gives  information  on  the  parameters  set  up  for  the  file  transfer  and  on 
the  progress  of  the  transfer. 

33  Implementation  Details 

This  section  presents  a  more  detailed  look  at  the  structure  of  TRANSFER.  Each 
major  segment  of  the  program  is  shown  in  a  flowchart  indicating  which  operations  are 
carried  out  in  the  main  program,  and  which  of  them  are  carried  out  in  subroutines. 


>  Figure  5  shows  the  file  selection  segment  of  the  program.  Wildcarding  and 

J  switches  are  implemented  in  this  segment.  The  user  is  first  prompted  for  a  filename,  and 

>  then  the  name  is  processed  to  check  for  switches.  If  switches  are  present,  they  are  decoded 

\  and  the  proper  flags  set  for  later  processing.  If  a  wildcard  character  appears,  a  flag  is  set 

for  tape  processing,  or  a  check  of  the  appropriate  directory  is  performed  and  all  matching 
names  extracted  for  disk  file  processing. 

\  Once  the  input  name  has  been  processed,  the  first  file  which  matches  all  criteria  is 

I  found  and  opened  (not  necessarily  an  easy  task  when  using  magnetic  tape).  If  the  file 

I  opening  was  successful,  the  DREA  header  is  read  into  a  buffer  for  use  in  setting  up  the 

v  transfer  parameters.  System  calls  are  used  for  all  of  the  magnetic  tape  operations  in  the 

£  interests  of  speed.  Disk  operations  are  done  using  a  set  of  subroutines  obtained  from  the 

>!  Naval  Research  Laboratory.  These  allow  asynchronous  operations  (ie.  computations  can 

•  carry  on  while  data  are  being  read  from  the  input  file)  and  are  written  in  VAX-MACRO,  so 

they  offer  a  speed  advantage  over  pure  FORTRAN  calls. 


When  a  file  has  been  properly  opened,  the  user  is  prompted  for  an  output  filename. 
A  default  (described  in  the  preceding  section)  is  presented,  but  this  can  be  changed  to 
anything  the  user  wants.  If  wildcard  files  are  used,  the  user  is  given  the  option  of 
specifying  transfer  parameters  separately  for  each  file  or  of  using  the  transfer  parameters  set 
for  the  first  file  for  all  of  the  others.  Output  format  (VAX  or  READRT)  is  then  chosen,  and 
the  output  file  is  opened. 

The  flowchart  shown  in  Figure  6  gives  more  detail  of  the  structure  in  the  processing 
and  data  output  sections  of  TRANSFER.  The  information  contained  in  the  DREA  header  is 
used  heavily  in  this  segment  of  the  program.  Parameters  such  as  record  size,  the 
multiplexing  type  and  the  number  of  channels  are  used  to  determine  the  the  number  of 
blocks  which  must  be  skipped  to  get  to  the  desired  point  in  the  input  file.  TRANSFER 
works  with  "frames"  of  data.  A  "frame"  is  defined  as  the  smallest  amount  of  data  which 
can  be  read  from  the  input  file  which  gives  a  full  set  of  input  data  for  each  channel  in  the 
file.  A  full  set  of  data  can  vary  from  a  single  time  sample  in  a  ".DAT"  file  to  a  number  of 
points  defined  by  the  FFT  length  used  for  processing  in  a  "  FTR"  or  ".PWR"  file.  This 
length  is  also  determined  by  the  type  of  multiplexing  used  in  the  file  (see  Figures  1,  2  and 
3).  As  well,  a  "frame"  defines  the  minimum  time  unit  which  can  be  accessed  by 
TRANSFER. 
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Determining  which  channels  are  to  be  processed  is  simple  enough  if  a  single  file  is 
being  used  as  input.  If  a  wildcard  file  set  has  been  chosen,  the  process  becomes  more 
difficult.  Channels  can  be  selected  by  "number",  that  is  by  their  standing  in  ue  order  in 
which  they  were  put  into  the  file  (first,  second,  etc.  in  the  multiplexing  hierarchy),  or  by 
the  hydrophone  channel  which  they  represent  (included  in  the  byte  label  of  the  DREA 
header).  For  example,  the  "first"  channel  recorded  in  a  ".DAT"  file  may  actually 
correspond  to  the  time  series  for  hydrophone  channel  #20.  The  user  has  the  option  of 
using  the  same  parameters  for  each  input  file  in  a  wildcard  set,  and  so  must  decide  whether 
to  key  on  channel  numbers  or  hydrophone  numbers.  While  such  a  set  is  being  processed, 
it  is  possible  that  channels  will  occur  in  a  different  order  in  some  of  the  files.  The  headers 
of  the  output  files  will  reflect  this,  but  it  is  often  better  to  key  on  hydrophone  numbers  in 
such  situations.  If  a  chosen  hydrophone  does  not  exist  in  one  of  the  input  files,  the  user  is 
notified  dunng  execunon  and  allowed  to  change  the  selection. 

Once  channels  have  been  selected,  the  DREA  header  is  modified  accordingly  and 
written  into  the  output  file.  At  this  point  the  data  transfer  can  begin.  Enough  blocks  arc- 
read  from  the  input  file  to  give  at  least  one  "frame"  of  data.  If  some  channels  are  not  being 
transferred,  these  are  then  removed  from  the  input  buffer.  It  is  most  efficient  to  write  large 
segments  of  data  to  disk,  so  input  blocks  are  processed  until  a  relatively  large  output  buffer 
is  filled  or  the  input  file  is  finished.  The  collected  (perhaps  reduced)  data  are  then  written  to 
the  output  disk  file  using  the  DBIO  routines  from  NRL  (for  a  VAX  format  file)  or 
FORTRAN  wnte  statements  'for  a  READRT  file). 

Once  a  file  is  completed,  input  and  output  files  are  closed  before  proceeding.  If 
w  ildcard  files  have  been  selected,  a  check  is  made  for  further  matches  and  if  any  are  found, 
the  next  matching  file  is  processed.  (The  user  may  or  may  not  be  prompted  for  transfer 
parameters  for  files  after  the  first,  depending  on  the  option  chosen.)  If  wildcard  files  are 
not  being  processed,  or  no  more  wildcard  matches  are  found,  the  user  is  asked  if  more  files 
are  to  be  processed.  At  this  point,  the  program  can  be  exited  or  a  new  run  started. 


4.  HOW  TO  USE  THE  PROGRAM 

A  sample  terminal  session  will  be  presented  in  this  section.  "Snapshots"  of  the 
user's  terminal  screen  w  ill  be  used  to  illustrate  the  user  interface  and  provide  information 
on  running  the  program. 

The  executable  version  of  TRANSFER  is  located  in  the  directors' 
SAS  IFARRELL.TRANSFER]  (soon  to  be  moved  to  DREAPACS:[TRANSFER]).  To 
stan  the  program,  the  user  must  type  RUN  SAS:[FARRELL.TRANSFER]TRANSFER. 

The  following  scenano  is  presented  as  an  example.  A  user  has  a  mag-tape 
(produced  by  a  PDP-1 1  data  recording  program)  containing  time  series  data  which  are  to  be 
transferred  to  the  VAX.  He  believes  the  filename  to  be  ZZZ003.DAT,  but  is  unsure  of  the 
extension.  He  does,  however,  know  that  the  file  begins  at  10:33:41  and  that  data  from  32 
channels  are  recorded  in  the  file.  Ten  seconds  of  data  beginning  at  10:33:50  and  4 
channels  out  of  32  in  the  data  file,  namely  2.  II,  18  and  23  are  to  be  transferred. 
TR  \NSFF.R  is  started  using  the  command  mentioned  in  the  preceding  paragraph,  and  then 
an  input  session  begins. 

In  the  following  figures,  a  header  categorizing  the  parameters  to  be  input  appears  at 
the  top  of  each  input  screen  (boldface  type).  After  every  prompt,  the  default  value  is 
presented  in  brackets  To  accept  the  default  value,  the  user  need  only  enter  a  carnage 
return  Help  can  be  obtained  after  any  prompt  by  typing  a  followed  bv  a  carnage 


return.  User-entered  responses  appear  in  boldface  type  (after  a  prompt,  so  there  should  be 
no  confusion  with  the  screen  headers).  Some  prompts  appear  only  under  special 
circumstances  (when  a  wildcard  file  input  set  is  being  used,  for  example),  and  these  are 
shown  in  italics  when  they  would  not  otherwise  appear  in  the  example  being  presented. 

The  first  screen  which  appears  once  TRANSFER  begins  execution  is  shown  in 
Figure  7.  Here  the  user  sets  up  the  input  and  output  file  specifications. 


File  Setup 

Enter  name  for  the  input  datafile.  (MSB0:Q38334.DAT?)  ? 

File  name:  UNIT  or  STRUCTURE:NAME. EXTENSION 

Enter  name  for  the  input  datafile.  (MSB0:Q38334.DAT  ?)  ?  MSAO:DUMMYDAT/V^) 

9cMOUNT-I-WRITELOCK,  volume  is  write  locked 

7cMOUNT-I-MOUNTED,  RT11A  mounted  on  _MSA0: 

Found  file:  ZZZ003.DAT 

Want  to  pass  over  this  file?(-l=yA>=only/l=no/2=rewind)  (-1?)  ?  1^) 

Default  output  file  extension  isTFR-  Change  it?  (N?)  ? 

Enter  the  new  default  file  extensi  on.  (TFR? )  ? 

Re-specify  output  parameters  for  e  ach  input  file?  (N?)  ? 

Enter  a  name  for  the  Disk  output  file.  (ZZZ003.TFR?) 

Should  the  output  file  be  Vax  (V)  or  READRT  (R)  format?  (V?)  ?  V 

NOTE:  signifies  a  carriage  return. 

Figure  7:  First  Input  Screen  of  TRANSFER 

To  demonstrate  the  use  of  the  "help"  feature,  a  "?"  was  entered  in  response  to  the 
prompt  asking  for  an  input  filename.  The  help  field  shows  that  the  name  format  is 
UNIT:NAME. EXTENSION,  and  the  prompt  is  then  repeated.  The  unit  or  directory  must 
be  the  first  entry  in  the  file  specification  -  if  no  unit  or  directory  is  specified,  the  current 
director)’  on  disk  is  used  as  a  default..  In  the  example,  the  tape  containing  the  data  is 
mounted  on  unit  MSAO:.  The  user  doesn't  know  the  name  of  the  file  to  be  transferred,  so  a 
dummy  name  has  been  entered  along  with  the  /V  switch  (this  means  that  the  program  will 
ask  before  skipping  any  file).  Messages  from  the  system  following  the  filename  entry 
show  that  the  tape  has  been  successfully  mounted  on  the  requested  unit  and  that  the  tape  is 
w  nte  locked  (ie  no  data  can  be  written  to  it  -  the  safest  policy  with  data  tapes). 

The  /V  switch  causes  the  name  of  the  first  file  found  (ZZZ003.DAT)  to  be 
displayed.  The  user  then  has  several  options.  If  a  -1  is  entered,  the  file  will  be  skipped 
and  the  search  for  DUMMY.DAT  will  continue.  Entering  0  will  cause  ZZZ003.DAT  to  be 
skipped  and  the  next  file  found  to  be  processed.  Entering  1  will  cause  ZZZ003.DAT  to  be 
chosen  as  the  input  file,  and  entering  2  will  rewind  the  tape  and  begin  searching  from  the 
start  of  the  tape  again.  The  default  is  to  skip  the  file,  but  in  this  case,  the  user  has  chosen  to 
process  it  by  entering  a  1. 


Next,  a  name  for  the  output  disk  file  must  be  specified.  The  default  filename  is  the 
same  as  the  input  filename  with  the  extension  TFR.  and  in  this  case,  the  user  has  chosen  to 
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go  with  the  default  by  entering  only  a  carriage  return.  (If  wildcard  input  filespecs  had  been 
used,  the  user  would  have  been  prompted  to  determine  if  the  .TFR  file  extension  should  be 
used  for  all  output  files.  At  that  point,  a  new  extension  could  have  been  specified  which  1 

would  appear  on  all  output  files  from  the  wildcard  set.  The  option  to  re-specify  output 
parameters  for  each  input  file  would  also  have  been  presented  if  a  wildcard  set  had  been 
chosen.)  Finally,  the  output  file  type  must  be  chosen.  Here,  the  user  has  chosen  VAX 
format  (described  earlier)  and  entered  a  V  followed  by  a  carriage  return  (  a  carriage  return 
would  have  been  enough  since  V  is  the  default,  but  entering  V  does  no  harm). 

The  next  screen  to  appear  deals  with  selecting  a  file  segment  to  be  transferred  and  is 
shown  in  Figure  8.  At  the  top  of  the  screen,  the  start  time  of  the  chosen  file  is  displayed. 

The  user  is  then  given  the  option  to  start  the  transfer  at  a  later  time  in  the  file  (default  is  to 
stan  at  the  beginning).  Here,  the  user  asks  to  stan  the  transfer  at  10:33:50.  TRANSFER 
skips  into  the  file  to  the  desired  point  (or  the  nearest  accessible  time  greater  than  that 
requested  -  constrained  by  the  input  file  format).  The  user  is  notified  of  the  stan  time 
which  the  program  is  actually  using  and  of  the  number  of  blocks  being  skipped.  The  next 
prompt  allows  the  user  to  set  the  amount  of  data  to  be  transferred.  The  user  has  three 
options;  n  (where  n  is  some  number)  will  transfer  n  512  byte  blocks  (to  the  nearest 
"frame")  into  the  output  file;  -1  will  transfer  all  data  from  the  specified  starting  position 
(time)  to  the  end  of  the  input  file  into  the  output  file;  and  -2  will  allow  the  user  to  specify  a 
time  interval  for  transfer.  In  this  example,  the  user  has  chosen  to  transfer  by  time  interval. 

Time  is  specified  in  the  HH:MM:SS  format,  and  in  the  example,  the  user  has  asked  for  10 
seconds  worth  of  data  to  he  put  into  the  output  file.  (If  a  wildcard  file  set  had  been 
specified  for  input,  the  user  would,  at  this  time,  be  asked  to  decide  whether  to  transfer  the 
same  channel  set  SSL  hydrophone  set  for  each  input-output  file  pair.  This  choice  doesn't 
deal  with  "data  segment  specification",  but  was  placed  on  this  screen  due  to  program 
structure  constraints.  The  prompt  appears  in  italics  as  wildcards  are  not  being  used  here.) 


Data  Segment  Specification 
This  file  starts  at  10:33:41 

Do  voc  want  to  begin  processing  at  some  other  time  (N?)  ?  Y  ■+) 

Enter  the  time  at  which  you  wish  to  start  ( 10:33:41?)  ?  10:33:50 

Actual  stan  time  will  be  10:33:50 
144  physical  blocks  will  be  skipped. 

Enter  n  to  X-fer  n  blocks.- 1  for  all. -2  to  specify  time  (-1?)  ?  -2 
Enter  the  length  of  ume  of  the  transfer  (00:01 :00r>)  0  00:00:10 

Select  tlu  same  H  P  <)'  i  or  channels  i.\i  from  each  file  (Y?  i  ? 


Figure  8:  Second  Input  Screen  for  TRANSFER 

Once  a  file  segment  has  been  defined,  a  channel  set-up  screen  appears  (shown  in 
F  igure  di  The  total  number  of  channels  in  the  input  file  (including  non-acoustic  channels) 
is  displayed,  followed  by  a  list  of  the  acoustic  channels  and  their  corresponding 
hydrophone  numbers  (from  the  byte  pan  of  the  DRF.A  headeri  From  these,  the  user 
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selects  which  channels  are  to  be  transferred  into  the  output  file.  The  number  of  channels  to 
be  transferred  is  specified  first  (-1  will  transfer  all  channels  -  including  non-acoustic  ones  - 
to  the  output  file),  and  then  the  channels  (not  hydrophone  numbers)  are  selected.  The 
selected  channels  can  be  separated  by  spaces  or  commas.  This  concludes  the  input  session 
as  TRANSFER  has  all  the  parameters  needed  for  execution. 


Channel  Setup 


There  art  32  channels  in  the  input  file. 


The  following  are  acoustic  channels: 

1,  2,  3,  4,  5,  6,  7,  8,  9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25 
27,28,29,30,31,32 


,26, 


The  corresponding  H/P  numbers  are: 

1,  2,  3,  4,  5,  6,  7,  8,  9,10,11,12,13,14,15,16,17,18,19,20.21,22,23,24,25,26, 
27,28,29,30,31,32 

How  many  channels  do  you  want  to  process?  (-1  for  all)  (-1  ?)  ?4 
Enter  the  channels  you  wish  to  study  :  2  11  18  23 


Figure  9:  Third  Input  Screen  for  TRANSFER 

Once  data  transfer  has  begun,  program  status  is  continuously  updated  on  the 
terminal  screen.  The  status  screen  (shown  in  Figure  10)  can  be  divided  into  three 
segments.  (Note  that  parameters  which  depend  on  choices  the  user  makes  during  the  input 
phase  of  the  program  appear  in  boldface  type  in  the  Figure.)  The  first  of  these  segments 
displays  the  current  time  and  date  and  is  updated  as  the  program  executes. 

The  second  segment  is  static  and  gives  information  on  the  input  and  output  files. 
Pan  of  the  ASCII  label  is  shown  to  give  some  indication  of  the  origins  of  the  file.  Below 
this  label,  the  names  and  stan  times  for  the  input  and  output  files  are  displayed.  The 
number  of  frames  requested  for  transfer  to  the  output  file  is  also  shown.  If  all  blocks  from 
the  specified  stan  time  to  the  end  of  the  input  file  are  to  be  transferred,  the  number  of 
frames  to  be  transferred  is  shown  as  'TO  EOF  which  stands  for  Jo  End  Of  £ile  (Note  that 
this  is  shown  in  the  sample  screen  of  Figure  10  -  which  came  from  a  different  run  of 
TRANSFER  with  the  same  input  file).  The  total  number  (including  non-acoustic)  of 
channels  in  the  input  file  and  the  number  being  transferred  are  shown,  along  with  their 
channel  numbers  and  hydrophone  numbers  (acoustic  channels). 

The  third  segment  is  dynamic  and  updates  as  the  program  executes.  The  number  of 
frames  processed  is  shown  along  with  the  average  time  taken  to  son  and  transfer  each 
frame  to  the  output  file.  In  the  sample  screen  shown,  100  blocks  have  been  processed,  and 
the  average  time  (clock  time)  per  frame  overall  has  been  .02  seconds.  A  running  count  of 
parity  errors  during  tape  reads,  is  given  below  the  "average  time  per  frame"  block  if  any 
errors  occur. 

When  all  data  have  been  transferred  to  the  output  file,  the  user  can  request  that 
another  file  be  transferred.  In  this  case,  the  program  starts  again  with  the  first  input  screen. 
If  a  wildcard  set  was  specified  as  input  and  the  same  transfer  parameters  were  to  be  used 
for  each  file,  the  next  transfer  will  proceed  automatically.  If  transfer  parameters  were  to  be 
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specified  for  each  file  in  a  wildcard  set,  the  user  goes  through  all  input  screens  except  the 
first  for  each  file.  If  a  tape  file  was  used  as  input,  the  tape  will  be  dismounted  and  rewound 
automatically  when  execution  of  the  program  finishes. 
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FILE  TRANSFER  STATISTICS 

4-DEC-1986  11:50:00 

ASCII  Label:  32  CHANNELS  REAL 

Input  File:  MSA0:ZZZ003J)AT 

Output  File:  ZZZ003.TFR 

Starts  at:  10:33:41 

Starts  at:  10:33:50 

Number  of  frames  requested:  TO  EOF 

Channel  usage:  4  chosen  out  of  32  total 

Channels:  2,11,18,23 

Acoustic  Channels:  2,11,18,23 

Number  of  Blocks  processed:  100 

Average  time  per  frame:  .02  sec 

Figure  10:  Program  Status  Screen 


Appendix  A  gives  a  list  of  all  prompts  which  can  appear  during  transfer  execution, 
along  with  the  'help'  string  for  each.  As  well,  full  descriptions  of  prompts  which  have  not 
been  discussed  fully  in  the  text  are  given.  Messages  which  appear  on  the  screen  when  an 
event  requiring  user  attention  occurs  during  TRANSFER  execution  are  also  listed  and 
explained  in  the  Appendix. 


For  completeness,  a  listing  of  TRANSFER  and  its  major  subroutines  is  included  in 
Appendix  B.  Some  of  the  VMS  system  calls  are  unavoidably  confusing;  however  the 
programs  are  fully  commented,  so  no  further  description  will  be  given  in  the  main  body  of 
this  document. 


5  FI  TI  RE  DEVELOPMENTS 


Currently,  output  files  cannot  be  written  to  magnetic  tape.  This  feature  would  be 
useful,  since  no  tool  for  reducing  the  size  of  data  files  exists  on  the  PDP-lls  (PDP-11 
user'  have  the  same  need  for  this  functionality  as  VAX  users).  Squeezing  several  ".DAT" 
fiies  onto  a  single  mag-tape  would  be  another  use  for  a  tape-output  feature  (no  tool  exists 
for  do.ng  this  on  the  VAXen  at  present).  Such  a  modification  requires  a  module  for  PDP- 
1 1  format  tape  output  (currently  being  developed  by  G.  Heard  at  DREA).  Once  the  module 
is  obtained,  the  upgrade  should  take  relatively  little  time. 


Transfer  of  other  file  formats  would  also  be  desirable  (for  example  READRT  to 
VAX  conversion),  and  will  be  implemented  if  enough  user  interest  is  shown.  A  generic 
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transfer  option  (block-for-block  copy  of  an  input  mag-tape  file  to  VAX  disk)  is  being 
developed  and  should  prove  useful  for  non-DREA  generated  tapes. 


6  CONCLUSIONS  AND  ACKNOWLEDGEMENTS 

This  note  has  described  a  robust  file  transfer  program  for  the  VAX  computer.  The 
program  allows  transfer  of  PDP-11  format  magnetic  tape  files  (and  VAX  disk  files)  to 
VAX  disk  files  which  are  formatted  properly  for  use  by  the  Surveillance  Acoustics  section 
suite  of  signal  processing  and  display  programs.  Data  manipulation  tools  which  allow 
selected  channels  and  data  segments  to  be  transferred  are  available  within  the  program, 
making  it  a  versatile  tool  for  pre-analysis  data  preparation. 

Thanks  to  Ed  Chaulk,  Vance  Crowe,  Phil  Staal  and  many  others  who  made  useful 
suggestions  during  TRANSFER'S  creation.  Almost  all  of  their  suggestions  have  been 
incorporated  in  the  working  program.  Asynchronous  disk  input/output  routines  written  by 
J.  Padgett  of  NRL  and  obtained  with  the  help  of  Art  Collier  at  DREA  proved  to  be 
extremely  useful,  increasing  the  speed  of  TRANSFER  significantly.  Subroutines  written 
by  Vance  Crowe,  Laurie  Bunch  and  Doug  Peters  were  also  used,  making  my  programming 
job  a  great  deal  easier.  Figures  1,  2  and  3  were  produced  by  Phil  Staal.  Thanks  should 
also  go  to  Bruce  Skinner  for  his  aid  with  VAX  VMS  system  calls. 


References 


Caldwell  D.A.,  "A  Standard  for  the  DREA  Data  Descriptor  Block",  D.R.E.A. 
Technical  Communication  87/302  ,  Dartmouth,  Nova  Scotia,  December  1986. 

Farrell  J.B.,  "SAFTR  -  A  Program  for  Producing  Fourier  Coefficient  Files  on  the 
VAX  Computer",  D.RE.A.  Technical  Communication  87/  (DRAFT) ,  Dartmouth, 
Nova  Scotia,  February  1987. 

Hurdle  B.B.,  "Private  Communication",  Naval  Research  Laboratory,  Washington, 
D.C.,  February,  1986. 


20 


APPENDIX  A 


A  -  DETAILS  OF  USER  PROMPTS  -  INTRODUCTION 

Appendix  A1  lists  all  of  the  prompts  which  can  appear  on  the  user's  screen  during 
TRANSFER  execution  and  gives  a  brief  description  of  each.  The  help  available  for  each 
prompt  is  also  listed.  Appendix  A2  lists  messages  which  appear  on  the  screen  when  an 
event  occurs  during  TRANSFER  execution  which  requires  attention  or  should  be  noted  by 
the  user. 

A1  -  Alphabetical  Listing  of  Prompts 

PROMPT:  Default  output  file  extension  is  .TFR  -  change  it? 

HELP:  All  output  files  will  have  the  specified  extension. 

DETAILS:  Allows  the  user  to  change  the  default  file  extension  which  will  be  used 

when  processing  wildcard  input  file  sets.  If  parameters  are  not  being 
changed  for  each  file  in  the  set,  output  filenames  will  be  the  same  as  the 
input  filenames,  but  with  the  default  extension  substituted  for  the  original. 

PROMPT:  Do  you  want  to  begin  processing  at  some  other  time? 

HELP:  Default  is  to  stan  at  the  time  shown. 

DETAILS:  Appears  when  the  stan  time  of  the  input  file  is  displayed.  The  user  can  stan 

the  output  file  at  the  same  time,  or  modify  the  stan  time  for  output  by 
responding  "Y"  to  this  prompt. 

PROMPT:  Do  you  want  to  process  another  file? 

HELP:  Default  is  to  exit  the  program. 

DETAILS:  Allows  processing  to  continue  when  the  current  input  file  (or  set)  has  been 

completed. 

PROMPT:  Enter  a  name  for  the  disk  output  file. 

HELP:  Default  will  be  the  same  as  the  input  name  with  the  extension  .TFR. 

DETAILS:  Defines  a  file  where  the  output  will  be  dumped. 

PROMPT:  Enter  n  to  X-fer  n  blocks,- 1  for  all, -2  to  specify  time. 

HELP:  -2  will  let  you  enter  a  time  interval  for  the  transfer. 

DETAILS:  Defines  the  segment  of  data  to  be  transferred.  Data  will  stan  at  the  specified 

stan  time  and  have  an  extent  specified  by  the  response  to  this  prompt. 

PROMPT:  Enter  name  for  the  input  data  file. 

HELP:  Filename:  UNIT  or  STRUCTUREtNAME. EXTENSION 

DETAILS:  Defines  a  file  (can  be  a  wildcard  file  set  and  have  switches  included  in  the 

filename)  where  the  input  data  are  to  be  found. 

PROMPT:  Enter  the  channels  you  wish  to  study. 

HELP:  No  help  appears  for  this  prompt. 

DETAILS:  Allows  input  of  a  vector  of  channel  #s  which  are  to  be  transferred.  This 

prompt  appears  after  the  number  of  channels  to  transfer  has  been  established 
and  a  list  of  the  available  channels  has  been  displayed. 

PROMPT:  Enter  the  length  of  time  of  the  transfer. 

HELP:  Format  is  HH:MM:SS. 

DETAILS:  The  user  must  enter  the  time  extent  of  the  dita  segment  to  be  transferred 

from  the  input  to  the  output  file. 
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Enter  the  new  default  file  extension. 

Typically  three  letters  long  (leave  out  the 

Allows  a  new  default  file  extension  to  be  set  Appears  after  the  prompt 
allowing  the  user  to  decide  whether  or  not  to  keep  the  .TFR  extension. 

Enter  the  time  at  which  you  wish  to  start 
Format  is  HH:MM:SS. 

Set  a  new  start  time  for  the  output  file.  The  prompt  appears  after  the  user 
requests  a  transfer  start  time  other  than  the  start  time  of  the  input  file. 

How  many  channels  do  you  want  to  process?  (-1  for  all) 

Enter  the  number  of  channels  to  process. 

Lets  the  user  select  from  the  channels  available  in  the  input  file. 

Proceed  using  the  subset  of  requested  phones  found? 

Re-specify  H/P  or  skip  this  file  if  the  reply  is  N. 

Prompt  appears  when  the  required  H/P  set  is  not  found  in  a  wildcard  file. 
Execution  can  continue  with  the  subset  of  H/P  found  in  the  file. 
Alternatively,  the  user  can  modify  the  H/P  set  or  skip  the  file  and  proceed  to 
the  next  one  in  the  set. 

Re-specify  output  parameters  for  each  input  file? 

Else  use  the  default  output  file  name  &  specs. 

An  option  used  for  wildcard  input  file  sets.  If  the  user  replies  "N"  the 
parameters  entered  for  the  first  file  in  the  set  will  be  used  for  all  of  the 
others,  and  the  output  filenames  will  be  the  input  names  with  the  default 
extension.  If  the  user  replies  "Y"  all  prompts  will  appear  for  each  input  file 
in  the  set. 

Select  the  same  H/P  (Y)  or  channels  (N)  from  each  file? 

Default  will  select  the  same  H/P  from  each  input  file. 

Used  in  conjunction  with  wildcard  file  sets  when  transfer  parameters  are  not 
being  changed  for  each  file  in  the  set.  Either  H/P  or  channels  will  be  kept 
the  same  for  each  output  file. 

Should  the  output  be  VAX  (V)  or  READRT  (R)  format? 

(V)  format  compatible  with  DISPVAX,  SASPEC,  etc. 

Sets  the  format  of  the  output  data  file. 

Skip  to  the  next  file  in  the  input  set? 

Otherwise  use  this  file  with  reduced  #  of  chans. 

Prompt  appears  (when  channels  are  being  kept  the  same  in  a  wildcard  set) 
when  channels  differ  from  those  expected  in  the  file  set  The  user  can 
proceed  with  a  transfer  of  the  reduced  number  of  channels,  or  skip  to  the 
next  file. 

Type  1  to  take  closest  record  start,  2  to  re-specify  time. 

Closest  may  be  earlier  or  later  than  the  chosen  time. 

Appears  when  the  user  has  requested  the  transfer  of  a  number  of  blocks 
which  doesn't  result  in  an  integral  number  of  data  frames  in  the  output  file. 
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PROMPT:  Want  to  pass  over  this  file?(-  l*y/0=only/l=no/2=rewmd  i 

HELP:  - 1  =pass  file&look  at  next/0=correct/l  =use  file/2=rewmd. 

DETAILS:  Used  in  conjunction  with  the  /V  switch  for  skipping  files  on  mag-tape.  The 

user  can:  (- 1 )  continue  searching  for  an  exact  match  for  the  input  filename. 
(0)  skip  to  the  next  file  and  process  it,  (1)  process  the  file  found  or  (2) 
rewind  the  tape  and  then  continue  the  search. 


A2  -  User  Messages  From  TRANSFER  (Grouped  By  Subroutine! 


A2.1  Messages  from  TRANSFER  main  program 

MESSAGE:  Stop  ***  fatal  --  input  file  is  empty. 

SOURCE:  TRANSFER  main  program. 

CAUSE:  Appears  when  the  user  tries  to  open  an  input  disk  file  which  contains  no 

data. 


MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 


Dismounting  Tape. 

TRANSFER  main  program. 

Tape  input  was  completed  or  end  of  tape  was  reached. 

A2.2  Messages  from  subroutine  Check 

Error  in  system  call  at  position  ***. 

Subroutine  CHECK. 

An  error  condition  resulted  during  a  VMS  system  call. 

A23  Messages  from  subroutine  MTFILE 

Stop  ***  All  files  processed. 

Subroutine  MTFILE. 

All  files  from  a  mag-tape  /START  /STOP  set  have  been  processed. 

End  of  tape  encountered  -  rewinding. 

Subroutine  MTFILE. 

Reached  the  end  of  a  mag-tape  being  used  for  input. 

Skipping  to  next  file. 

Subroutine  MTFILE. 

A  file  with  a  name  different  than  the  input  filename  was  found  and  is  being 
skipped  (either  because  the  user  requested  it  or  because  the  /V  switch  was 
not  used). 

File  ****  found. 

Subroutine  MTFILE. 

Found  a  file  which  matches  the  input  filespec. 

Found  file  ****. 

Subroutine  MTFILE. 

Appears  for  every  file  found  when  the  /V  switch  is  used  in  the  input 
filespec. 
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MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 
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Assume  you  want  to  skip  to  the  next  file. 

Subroutine  MTFTLE. 

Appears  when  the  user  responds  ambiguously  to  the  prompt  asking  whether 
or  not  to  skip  the  file  found. 

A2.4  Messages  from  subroutine  GETTAP 

Tape  is  already  mounted  -  assuming  you  did  it 
Subroutine  GETTAP. 

A  mount  request  has  been  issued  for  a  tape  unit  that  was  already  mounted. 

Error  -  device  probably  allocated  to  another  user. 

Subroutine  GETTAP. 

Failure  in  a  request  to  allocate  a  tape  unit  to  the  TRANSFER  job. 

Device  already  allocated  to  you. 

Subroutine  GETTAP. 

Routine  tried  to  allocate  a  tape  unit  already  allocated  to  the  user. 

A25  Messages  from  subroutine  Channd  Sdect 

Channel  Select  -  no  acoustic  data  found  in  file. 

Subroutine  ChanneLSelect. 

The  byte  label  of  the  DREA  header  indicated  that  there  were  no  acoustic 
channels  present  in  the  file. 

Working  on  file  ***. 

Subroutine  ChanneLSelect. 

Informs  the  user  of  the  current  file  being  processed  in  a  wildcard  file  set. 

The  H/P  available  differ  from  the  originals. 

Subroutine  ChanneLSelect. 

One  or  more  of  the  requested  hydrophone  channels  is  not  present  in  the 
wildcard  file  being  processed. 

You  are  keying  on  channels  rather  than  phones,  so  I  am  proceeding. 
Subroutine  Channel_select 

Hydrophone  numbers  have  changed  in  a  wildcard  file-set  but  the  user  is 
keying  on  channels  rather  than  phones,  so  this  serves  as  a  warning. 

You  asked  for  a  channel  not  found  in  the  input  file. 

Subroutine  ChanneLSelect. 

User  is  keying  on  channels,  and  one  (or  more)  of  the  required  ones  was  not 
found  in  the  current  wildcard  input  file. 

Proceeding  with  reduced  #  of  channels. 

Subroutine  ChanneLSelect. 

User  has  chosen  to  process  a  file,  even  though  it  doesn't  contain  all  of  the 
files  originally  asked  for. 

There  are  ##  channels  in  the  input  file. 

Subroutine  ChanneLSelect. 

Informs  the  user  of  the  number  of  channels  available 
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MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 


MESSAGE: 

SOURCE: 

CAUSE: 

MESSAGE: 

SOURCE: 

CAUSE: 
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The  following  arc  acoustic  channels. 

Subroutine  Channel_Select.  * 

Informs  the  user  of  the  acoustic  channels  in  the  input  file. 

The  corresponding  H/P  numbers  are: 

Subroutine  Channel_Select. 

Informs  the  user  of  hydrophone  numbers  corresponding  to  channel 
numbers. 

Channel  Select  -  sorry  no  default 
Subroutine  Channel_Select 

The  user  selected  no  channels  for  transfer  to  the  output  file. 

A2.6  Messages  from  subroutine  Skipper 

This  file  starts  at :  HH:MM:SS. 

Subroutine  Skipper. 

Informs  the  user  of  the  start  time  of  the  current  input  file. 

Actual  stan  time  will  be  HH:MM:SS. 

Subroutine  Skipper. 

Informs  the  user  of  the  actual  stan  time  of  the  output  data  file.  May  be 
different  from  the  requested  time  due  to  data  frame  size  in  the  input  file. 

Disk  start  block  will  be:  ###. 

Subroutine  Skipper. 

Informs  the  user  of  the  start  block  number  of  the  input  disk  transfer. 

###  physical  blocks  will  be  skipped. 

Subroutine  Skipper. 

Informs  the  user  of  the  blocks  to  be  skipped  on  tape. 

You  specified  a  zero-length  transfer  -  try  again. 

Subroutine  Skipper. 

User  specified  a  transfer  time  less  than  the  frame  time  of  the  input  file. 

WARNING  -  you  are  trying  to  transfer  0  blocks  -  try  again. 

Subroutine  Skipper. 

User  specified  a  number  of  blocks  less  than  the  number  of  blocks  in  a  data 
frame. 

Transfer  time  must  be  at  least:  HH:MM:SS. 

Subroutine  Skipper. 

Informs  the  user  of  the  minimum  time  which  can  be  specified  for  a  transfer. 

You  are  not  using  an  integral  number  of  records. 

Subroutine  Skipper. 

User  tried  to  transfer  data  which  did  not  fit  into  an  integral  number  of 
frames. 
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A2.7  Messages  from  subroutine  Disk_Wildcard 

MESSAGE:  No  more  files  match  the  input  spec. 

SOURCE:  Subroutine  Disk_Wildcard. 

CAUSE:  All  files  matching  a  wildcard  specification  have  been  found. 

A2.8  Messages  from  subroutine  Reader 

MESSAGE:  Error  reading  record! !  Error  count  =  ##. 

SOURCE:  Subroutine  Reader. 

CAUSE:  Routine  encountered  a  parity  error  on  tape  and  successfully  passed  it 

MESSAGE:  End  of  file  encountered. 

SOURCE:  Subroutine  Reader. 

CAUSE:  End  of  an  input  file  was  encountered  during  a  read. 

MESSAGE:  Saving  ###  frames  and  exiting. 

SOURCE:  Subroutine  Reader. 

CAUSE:  Informs  the  user  how  many  data  frames  are  being  stored  after  an  end  of 

file  was  encountered  during  a  read. 

MESSAGE:  Read  puts  us  past  EOF  -  blocks  to  read  ###. 

SOURCE:  Subroutine  Reader. 

CAUSE:  Disk  read  would  go  past  EOF. 

MESSAGE:  New  blocks  to  read  -  ###. 

SOURCE:  Subroutine  Reader. 

CAUSE:  Informs  user  of  the  number  of  blocks  which  will  be  read  from  a  disk  file 

when  a  full  read  would  go  past  the  EOF. 

MESSAGE:  Error  encountered  on  disk  read  -  saving  what  I  can  and  exiting. 
SOURCE:  Subroutine  Reader. 

CAUSE:  Informs  the  user  that  an  error  occurred  during  a  disk  read. 
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B  -  PROGRAM  LISTING  -  INTRODUCTION 

This  Appendix  contains  listings  of  the  TRANSFER  program  and  us  major 
subroutines.  The  casual  user  is  unlikely  to  Find  these  listings  useful,  however  those 
writing  programs  for  accessing  PDP-11  tape  Files  should  be  able  to  glean  some 
information.  The  programs  are  commented  fully,  and  most  operations  other  than  system 
calls  should  be  easy  to  interpret.  Print  size  has  been  reduced  to  conserve  space. 

B1  -  Listing  of  the  TRANSFER  Main  Program 


PROCRAW  NAME:  TRANSFER 


Written  by: 
Joseph  B.  Farrell 
DREA 

21  Jan.  1986 


Latest  revision:  1  Aug.  1986 

This  program  accepts  a  fl  lename,  opens  t  nat  file  cr.  ar.  RT-Il 
rr.ag  tape  (cr  VAX  disk  file  in  standard  'PDPllS*  format;  and  reads  ir.  t he 
DREA  standard  header.  The  user  car.  then  choose  some  or  all  of  the  input 
channels  tc  c>e  written  lntc  ar.  output  distt  file  (ir.  standard  ‘PDPllS’ 
format  cr  ir.  the  old  * READRT 1  format;.  The  tape  unit  does  not  need  to 
De  mounted  or  allocated  before  the  program  Is  run. 

The  program  was  written  tc  replace  anc  expand  or.  the  funct Iona.  Ity  of 
the  READRT  program  written  by  Ker,  Hahn  of  ASP. 


C  SUBROUTINE  CALLS: 

c 

I) 

ERASE  SCREEN- 

c 

luri  tn«  acrear.  or  tut  uaer’s  tanalna. 

c 

2) 

VTMESS 

Pott  a  message  or.  tna  user '  >  tenr.ina. 

c 

3) 

INPUTS 

Raaai  a  string  from  -jeer's  termina. 

c 

4! 

SYSSASSI3N 

Subroutine  to  assign  a  cnanna. . 

c 

5) 

LIBSSTOP 

Stops  program  execution  or.  arror. 

c 

6) 

SYSSQI  Or 

Queued  input-output  iron  a  channel. 

c 

7) 

CHECK 

enact  m  tn*  atatua  bc!!*r  a  .'tar  a  CTO 

c 

8) 

GETYN 

Gats  a  Yes  or  No  response  from  tna  uaar. 

c 

9) 

D  BO  PEN  (C) 

Opart  (craataa)  a  disk  Ilia. 

c 

C) 

DBSIZE 

Dataminaa  tna  tits  o!  a  dlak  file. 

c 

1) 

D  BREAD 

Raaaa  data  fro*  a  disk  flia. 

c 

2) 

DBWAIT 

Walt*  for  disk  operation  to  complete. 

c 

3) 

DBKUTE 

Turns  of  arror  reporting  from  DBKT  routlnas 

c 

4) 

SKIPPER 

Skips  into  a  fi.a  tc  a  apacifiao  time. 

c 

5) 

Channe._ select 

Cnoosts  cnan.na.s  to  procaaa. 

c 

6: 

Write  header  - 

Writes  header  to  the  output  file. 

c 

7) 

inputT 

Reads  ar.  Integer  from  the  user. 

c 

8) 

DBWR1TE 

- 

Writes  data  tc  a  cisic  file. 

c 

9) 

DBCLOSE 

- 

Closes  disk  files. 

C  20) 

SET  CURSOR 

• 

Moves  the  cursor  to  s  specified  .ocatlor.. 

♦  * 

PROGRAM  TRANSFER 


Parameter  and  Internal  variable  declarations. 
IMPLICIT  INTEGER* 4  <a-z) 

INCLUDE  • (Sdmtdef ) • 


PARAMETER  MSG-5 

PARAMETER  ( SSSENDOFTAPE- 1 §70 • XI 


BYTE 

BYTE 

BYTE 

BYTE 

BYTE 

BYTE 


biabel (128; 
s:ore_b.abel (128) 
EOT_search (1C) 
output_data (1C24C) 
raw_byte_data (5CC0C) 
raw  header  (512) 


CHARACTER 

CHARACTER 

CHARACTER 

CHARACTER 

CHARACTER 

CHARACTER 

CHARACTER 

CHARACTER 


a .aoe. *32C 

de!a«.it_extensior.*64 
char.  *6 
defaul t  *1 
lfl 1 e*  64 

newt lme*fi 
oidt lme*8 
of;ie*64 


Byte  part  of  DREA  header 
:Temp.  store  for  byte  *abe* 
: Checks  for  er.o  of  tape 
'Output  data 
Raw  input  data 
Raw  header  data 

ASCII  part  of  DREA  header 
■  Default  fl.e  extension 
1’Jaed  in  mag  tape  channel 
:i’ser  Input  defau.t 
•Input  filename 
Start  time  for  processing 
Start  time  of  Input  data 
'Output  filename 
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CHARACTER 

otypa 

'F.ag  for  output  f  1  .•  typa 

CHARACTER 

lApi  juri  i  1C : 

Cnacfcs  t*p*  uru 

character 

tip*  n**a '1C' 

Haaoar  filtniM 

CHARACTER 

wf i.#*6* 

WorKlrvc  fl  .*n«M 

INTEGER* I 

d.  ock  »_  f  ***,.  1  nput 

B . ock *  to  raaa  fro*  fli* 

INTEGER  *1 

pyta*  p»r_t’’  c* 

Dat*m nac  py  input  fi.a 

integer*; 

fi.a_fyp*~ 

Input  fliatypa  (.DAT.  ttc.i 

INTEGER*; 

nydropnona* (1 20 

Hyaropnona  maun 

integer*; 

: .  *o* .  (16, 

Intagar  part  of  DR LA  n«ad«: 

INTEGER*; 

input  cnanna. 

Input  cnanna.  for  tapa 

INTEGER*; 

r.  u»par_c!_cnar.na.  a 

#  channa.i  in  input  fi*a 

INTEGER*! 

ni«par~aort.ad 

NiAtxr  of  *I«Mnts  sort  ad 

iktege  ; 

procaas (120? 

Cnanna. »  tc  ba  procassac 

1  sTEG*.  •: 

Mi*ctor_Mi(  i!20i 

Tor  picking  bytas 

INTEGER*; 

•ift_f .ag 

f.ag  for  channa.  ramova- 

INTEGER*: 

■tart  fi.a 

Start  for  numoarac  tapa  fi.ai 

INTEGER** 

■top  ?:.* 

Stop  for  numoarac.  tapa  fi.a» 

INTEGER*. 

*.•»*.  ~lCSt  1  * 

Status  ouffar 

I NTECER  *  < 

c.ocr*  par_racera 

Oat  anti  nac  py  input  fi.a 

I  NT  E  GE  r*< 

CiOCki~t  c_ raaa 

#  block*  for  .nput  raac 

integer** 

cicc*»2#-  o’ writ* 

#  P.-cks  for  output  writ# 

INTEGER** 

Dyt •  i  :c_;m3 

#  oytas  raad  oaf  ora  aort 

integer** 

c.oar.na.  *”t  c  proc*** 

t  c na.nna.s  to  procaas 

I  NT  EGE  R  *  4 

ccnsa cut  Tva~cr.anna .  Dytat 

INTEGER *< 

3111 

F.ag  for  disk  fl.a  input 

INTEGER** 

d„»k  cnanna. _out 

’Output  cnanna. 

INTEGER** 

Irm  t:i* 

i  cl k a  or  racorot  par  fraaw 

INTEGER** 

frarr'aa^tc  raac 

input  f raMi  to  raac 

INTEGER** 

1  D.  « 

B.oek  for  disk  raada 

INTEGER** 

lar ror 

Error  fror  dak  raadt 

INTEGER** 

It  ima 

INTEGER** 

a:»*  input  _cnanr.a. 

Input  cnanna.  for  disk 

INTEGER  ** 

*  .rip ' 

incrawant  for  aort  mg  cnanna. 

INTEGER** 

f  rai*a»  *t  crae 

countar  for  framas  atorac 

INTEGER** 

*  rat 

cour.tar  for  fraina*  prouaaaac 

INTEGER** 

ma»«  : 

INTEGER** 

no. ocu  * 

NuJTiDar  of  d.ocij  :r.  input  fi. 

INTEGER** 

OC. OCX  * 

dafau.t  i  of  output  D.ocka 

INTEGER** 

of!  «r. 

uaaa  for  cnanna;  sorting 

INTEGER** 

*y*5a:  »jtk3- 

C: amount*  *  tapa 

INTEGER** 

tD.  OCK» 

INTEGER** 

tot*.  fc.OCk* 

ldgi ca: 

■  not  nar 

procaa*  a net nar  f:»a. 

logical 

ar>c_cf  tap# 

Enc  cf  tapa  indicator 

cog i cal 

•*t  r. 

F.ag  for  output  aatanaior 

log: cal 

:orr*c'  fi.a 

Fi.a  f.ag 

LOGICAL 

3*‘  y 

Laar  raaponaa  furct.or 

LOGICAL 

out  apa  c 

Ra-oc  output  fi.aapaca 

logila: 

»«-C 

Sup  Ilia  f.ag 

i,og:ca: 

var 1 fy 

*  a  r  i  f  y  fi.aa  tc  t>a  akippac 

log:  cal 

m i . dears 

*..acarc  .r  fi.anama 

RiAl 

c .  oc»  ; lx* 

R  £A- 

L.tir 

d.ocks  tc  raac 

sea: 

C.«Sw 

d.ock*  to  wrJta 

REAL 

Pi  kit 

tamporary  atoraga 

REAL 

f .aoa.  1 1 

F. oat  mg  part  of  DREA  naadar 

REAL 

fracticr.a.  p.ocki 

REA! 

frajna  fraction 

REAL 

frana  tima 

REAL 

nrac 

REAL 

racorca  par  raad 

RIAL 

sair.p.aa  par  input  c.oci 

Data  nr  mad  Dy  input  fi.a 

REAL 

t  '.m*  i  ntarva . 

1 v* .*- :f 

tapa  r.ama  1 1  r«^  naaoar  ri 

aq-i v« . *' c« 

.  . *  Da .  ra»  naaaa-  I 

•  q - 1  v *  .9’  :~r 

LC"  aaarcn.'apa  na:« 

f.aba.  .  ra*«  “aadarf33 

•  q- l v* . a- re 

t.aoa.  tmw  naadar<fc^ 

•  q  .  I  v*  .  *r 

*.ar>a.  ,  ra*  naadar  19: 

:ofl»c :  ic  *  1 

a t  1  ■  *.  1  c a  *rac.nrac 

rorro*  flTar 

wf : . a, cf i . « 

r  ojmo  -  o  ca  t  * 

o  jt  p  .  *-  a*  t  • 

;o imor  raw  r 

■to  pyt#  dat * 

aatarna. 

5  **:pf..a,.r.f  »  .  :5  i« 

C r ace  :  c .  i  5  r .  *  c 

■  *.  *  fc:  .it: 

a-  *  :  r  i  *  •  •  *sp  :  n<  .::A 

*. :  ma  -  . 

*  D - OCk t  - 

»:ir'  f:.*  - 
■top  ?!.#  -  9  999 
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default  extension  •  * TFR ' 


a r  t n«  screen  ana  writ*  4  n*aa*r  mimq*. 

call  *r*e*_scr**n (i, I » 
c4 . 1  s*t _cur *or (2, 1 J 

C4  *  *  vtmiss f • r*' , ' c* , ‘  Fii*  Setup  ’  VT100/20C 


Stt  4  name  for  th*  input  fli*. 

C4 . .  Inputs  (prompt  (1 ) ,  1  f lie,  h*lp  (l  i  ! 
C4..  StrSupCAS* liftl*  ifli*i 


D*t*nrin*  what  swicn*s  «*r*  us*c  ir.  the  input  film  specification, 
cs . .  switch***  :  f  i  i.*,  v*rl  fy,  start  _f  •;*,  st  op_fi  i*> 


A33  a  .DAT  extension.  if  non*  was  supp.  i*c  and  tn*  fi**r.am*  was  no 
sp*c.fl*c  as  a  »i.acarc, 

-i.acaro  •  .FALSE. 

1 f  ;  ;no*x  1: f; .n*. Cl wi.ocarc  -  .TRUE, 
if  1  md*x  if;.*.  ■  .  ’)  .*q.Cith*r. 

1  -  . *n2  ' i f !*•: 

1  f  ( 1  f  1 1  •  ■:  1  1  .  n*  .  *  •  *  >  t  n*r. 

.fl**  -  If  1 *•  { :  i*r.2  *  1  f  1  i*»  j  '  .  DAT  ’ 

•  nc  if 
«nc  1  f 

l!  (start  f : **. n*. T: wi .dcaro  -  .TRUE. 


Tree*  it  s*«  if  it’s  a  o.st  or  a  magtape  fii*. 

3is«  -  : 

.  f  :  nde*  i :  f . .  * .  .  eg .  Z  .  cr .  :  1  f : . •  ( I  ;  2 :  .  n* .  ‘MS  * . and . 

.  f  :.*”.:  2  .  n*  .  *«T  4  3.  s«  - 


if  a.s*  •  :  gc  f  n*  c: as  file  op*r  sectio',  ctn*rwis*  proc**a  wi 
«  rragtap*  op*'  ar.c  r*ac. 

.  f  •  3.  s«  .  n* .  1 :  t  r.*r 

man-  1  f  1  .#  1  :  5  : 

ta . .  mt f  .  .  •  input  cnarr*.,  raw  r«aot:,  if:.*,  wf:.*,:t:»*. 
v*r  :  f  >  .  start  f  1 . • ,  tt  op_  f  1  T* 

*.s* 


-■per  arc  :  *ac  fror  t  n*  a  is*  f..*  if  tap*  .5  not  tn*  input  medux. 


f 1 .*.  it im*. wf i .* 

-  if:.*'.,  .er.2  !  1  f  l  .*: 


cr.ar.n*.  ,  -  f  i  .* 
op*' 


e  •  •  m. .  n*  :  r*  jrx>* :  c  f  r  .  o c  *  s  : '  t  n*  i  np„t  f  1  .  *  . 

-4 .  .  dL»  .  j*  *c  1  s*  _  :  np*t  _ cr.ar.n* .  ,ns.c:u 
.  f  'ht.ocss  .eg.  C  S' ot  •••  fata.  --  input  fi,#  is  empty 


Reas  in*  n#ao*r  fro r  in*  input  fi.*. 

cs  •  -  dcread  1 31 *«  1  npur  cr.ann*. ,  raw  reader,  1 
ua..  aowa  1  ’  ■  a:s«_  .rput  cr.ann*. 


'  *  *  ry  •  .at1*  .'*c  'erptrary  s*  crag*  fir  .a**:  mod  .  f  .  cat  .  or. . 

3-  -  .  12* 

s'-  - r *  .  -  c.ao* 

•  nc  3 : 


«-  "•#  si'*#'  a '!  4  *m  .message  fnr  wi.dcajc  f:.*s  aft*r 


1'  .  . m*  g*  . .  i’ :  *. .1:1: 

ca . .  *ra s*  s'r «*r 

a  .  .  t*'  i'j.'5'  r  « ,  i 

a  •.".-* ss  ■  '  re  '  .  '  c  '  . 
*n  :  .  f 


•  '  '•  - s* r  '.c  arc*pt  arc  cr  c range  tre  defa-.t  cutpu*  fi  .erant. 

-f  .  m*  .  eg  .  .  .  and  .  .  acardi  t  r*- 

d*  f a_ . t  -  K  ‘ 


c  a'C .  c„t  spec  *.  '** 
r  1  .e  s«:  j; 


1  f  (*. :  .acard  t  n*r 

c 4 . .  cis*  -1  -dcara  1 

*  .se 

wfi.*  l:.*r.i  if:.* 

•  nc  :f 

ca  . .  demy  •  i  b 

ca..  doep*-  ■c:a«_  i-p-.t  _ 

. *  i *r : . n* . .  t  n*r 

t  >*p* •  .  Error  cr  f .  .* 
c  .  s*  _  :  r.p-t  era  .on*.  - 
. *r :  ■  1 


»r.c  1  f 


SPECIF 


APPENDIX  B 


APPENDIX  B 


axt  r.  -  gatyr 

lyafa-.t  output  ; . .  •  Axtar.sic-  .  *  tfr  -  cnanga  :t  , 

K.  .  output  fl.«*  *:..  nava  tna  spacif.a:!  AAtar.»;cr'. 
da ! a u *  t  : 

:  !  <axt  r.;  t  nan 
ca..  input* 

Er.tar  tna  na»  dafa..*  f:.t  axt  an tier  .  . 
aaf *u . :  ^axt  ansior., 

?ypTca.*>  thraa  .attar*  .cr.g  utivi  cut  th#  . 

and  1 f 

:  f  im.  .dcard;  mtr. 
dafa-.t  -  ’S' 
ojt*p«c  -  gat  yr 

fca-spacify  output  param#-.  an  fer  aacr  input  !:*•? 

L.sa  -se  "i*  oa fa-.'.  o-.put  !  : .  ar-aiaa  .  ‘  , 


r.a  daia-.t  output 
.  t*  •  r.*  Axtansic- 


fi.Ar.ariA  tc  da  t  na  sara  4* 
dafa-.t  axtans.o- 


input  f  .ar.a 


.  -  ;no»x <w! : . a. ' . 

:f'..aq.C.:  -  .ar.2(wf:.A 

c  !  i  .  a  •  wf  1  .a  iiu«  :  not*  i *f 1 . a .  '  -  1  ,  k  i  no«*  Ixfl**,  1  ‘  •  1  ; 

oafa-.t  axt ansi  or  •  1  :  .  ar.2  *aaf a».t  axtanalon:  ! 

:  f  •  it ina.aq. 1  .or  . out spac i  tnar 


*.  na  .*•: '  1  cncica  Jo:  tna  output  fl.anama. 


ca .  .  1  r.put  s  (prompt  *3 


rirt'.n  wnr.f*:  t na  output  f:.a  .* 


da  vax  c:  READ??  format  . 


ot  ypa  -  '  v  • 

ca. .  input* (prompt  (5.  , ctypa, na.p d. 
ca..  st rSupcasa  lotypa, ot ypa 


upa-  t  f.a  f..«  if  tna  raq^ASta::  forma*.  . 

IBC? ENC  .»  4  dis*  I  *•  C  routine  fror  NRl.  cr.g.r.a.y  written  for 
if.  mAcrc. 


. f  lotypa .•q. *V‘ >  tntr 
ca..  dtsnutad: 

ca.  .  dDopar.c  <di  tx  cnannt ._out ,  0!  1  .a1 

•  .it 

op#r.  (unit  -2,  fi  . a-of i  .A,  *t  At  us-‘ haw  •  .  form-  *  urformattad  ’  , 
*  ACCAAA- • saquant i a. * 

and  if 


"cr«  out  pA r Atn#t tri  using  3a t a  fror  tna  input  haadar. 


f 1 :  •: ypA- : 

:  :  '  (I.ADA.  r.  C  '  .  AQ  .  2  .  .  AP.3  .  <1  .ADA.  (11  !  .  AQ  .  1  !  >  f  i  i#_typ*-2 
;  f  '  '  i  *ada  .  1 1 C .  .  eq .  2  .Ana .  C  .ada.  i  1 1 ,  .  na .  i  ’  ;  f  i  .e~t ypa- 3 
t -r^>Ar_cf _chAnnA.s  -  ilaoa. d; 

D/'t  aspar  ~c.ock  -  i.AbA.d;  •  2 
c.oc«R_pAf_rAcora  -  : .ada. (2 <  ' i Iada 1  d ; 

*Arp.A*  par  ir.put_b.ocK  -  2  .  *  f  .  oAt  ( i  2  ada  id)* 
f . oAt  (iTaoaT  <fc;  T.umDAr  of  cnannalt) 
fract ion*.  D.ocks  -  floatloytA*  par  Diock; ■ 
float  ( r. umba r  _cf_c ha nna  Ts  •  i  Iada.  (6)  ) 
f  riict  ;r:-.A.  blocks  -  f ract 1 ona . _bl ocks  -  int  ( 
f  fact  lona  l__biocx*) 

:  f  f  f rac: :  ona.  _b.  ocu. i t . C. COOCI ) t  nan 
b.ocK*_p«r”*At  -  1 

A  .  AA 

b.ocK*_pAr_»At  -  int  ( (1 . /fract :ona._block») ♦().  5) 
an  c  1  f 

f ramA_f ract 1  or  -  2 . *f.oat (I.adaI  d ) ) •  (f .oat  (iiabal (6) ) • 
float  ( i i ada .  f.oat  (I.ada.  (8)  )  : 

. f  i  f 1 .e_typA. aq. 2) f rame_t lma  -  f*aDe.(B: 

Ififi.e  typA.aq. 3, f ramat lma  -  f*abA.(6) 

L.ccK_tImA  -  frama  fractlor,  •  frama  tima 
framA~*iiA  -  1 f lx  (7* oat  (i .aba.  (S )j •Tioat (1 .ada.  P) ) • 
float (1 1 ada 1 (6, ) /float (Dyt A*_pAr_blocM ) ; 

1  f  !  frAma_6l  z*  .  .t  .b.ock*_pAr_rAcord)  tnar. 
if  (diiK .aq. 1 .and. frama_si za. it . 1) than 
b.ocKS_tc_raacl  -  1 

A  1  SA 


blocK*_to_raaa  -  Diock*_par_rAcora 
and  if 


A  !•  A 

fcloc«*_tc_raaci  -  framA_*lzA 
and  1  f 


Ca..  routina  wf.ich  o«'»  tna  aasirad  starting  tima  and  staps  Into  tha 
fl.a  to  tnat  tlma.  Tha  routina  a. so  dAtanrinas  tna  total  numDar  of 
t*oc*e  to  da  raad  from  tha  Input  fiia. 


call  ik ippar ( input ^channel , samp. as_par_ input _b.ocK, f *aDal { 1 ) , 
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a  .  aoe  . ,  c.  oc*s_per_recorc,  c  .at  :m«  ,  newt  ime,  c:  ** ,  ibis, 

f  1  *e_type,  out  spec,  ;t  ;me,  tota._oioc*s,  o*  ocR_t  ime, 
f r ame _ t  :»• ,  fract:ona._c.oc*s,  numoerof ^channe. s j 


>*  a..o*  the  user  to  choose  the  numoer  cf  channels  to  put  into 
.e  output  file 


ca..  channe.  select  {D.aoe. ,  Nujr.i>er_c;_channeis ,  out  spec,  it  ime, 

Channe.  s  tc^process,  proceasThyaropnones,  Sift_flag,wfii< 

w;idcara7**ip_f i*« 

:f is*ip_fi.e.eQ. 1 > go  tc  1 C2 


id.iy  t.ne  header  tc  reflect  changes. 


: .aoe .  :  1  •  2 5€ 

:f  -l.aoe.  (8;  .eq.lither. 


i«t  sure  the  recorc  size  ir.  the  header  Is  ar  Integra,  number  cf 

...OCXS  . 


:.ar>e.  !2  -  aoe. •!' 'channe. s_t c  process  '  n  ujrx>er_of_cr.anne.  s 

;  f  ir.oc  ; .aDe.  (2:  ,  2*8.  . ne . C.  i  *aoe.  (2>  -256 
end  if 

i.aoe.  ■*  -  cr.anne: s_tc ^process 

ifii.aoe.  -;3  .  ie  .  C>  i  .aoe.  ;3  -  -1 


;rt  the  process  vector 

oc  *  -1 .  (Channels  tc_process-l ; 

ar  *  •  (K-l ,  ,  cHar.nels_tc_process 
. f  .process  ‘  . it .process (t :: t her 
hc.c  •  process  (x . 
process (x,  -  process:^ 
process  ( *  •  he : c 

enc  i  f 
enc  do 
er.c  ac 


at  ;  -  ,  channe .  s__:  o_proeess 

*  -  (« 'process  \*T'>  ~  1 
c.aoe.i.,  -  store&.aDe. ix' 
c.aoe. ( . *  1 >  •  store  b.aM.  U«'.i 
.  -  .  ♦  2 
•r.a  dc 

ac  *- f  i2*cr.anne:  s_tc_process  .  *  .j,12e 
D^aoe. { * :  •  C 
end  dc 


the  selector  mass 


vector . 


k  •  1 
:*  -  : 

ac  ;  •  1 ,  njnoer^cf _chahne.s 
.  !  i  ‘x  .  eg  .  process  (X1.  )  tner 
se.ect cr  masx  1  *.  •  *. 

«  -  *  *  T 
e .  se 

se.er.  ©r_rr.asx  <  4  i  •  C 

er.o  : ! 

'k  •  *  >.  •  1 

enc  dr 


:  _^ta  into  the  output  f:.e 
d .  '  V  i  t  her. 

.  t  eid.su  channel  out .  ;*>.  header,  1  j 
-a.t  (a; sk _ cKsnne . _out > 

2  « raw_ header  ( i j  , l -1 , 51 2 . 

*  '-.arne.  _eytes  -  :.aoe.»6  •i.aoe.'B 
•  *•  r-t :  ve_c.nar.ne.  _oyt  es  *  r.ur-Der  c*  channe.  s 
e  :  -  r . o c*s_t f _ reac  •  oytes  per_tIocK 
fT  reac  •  !.cit  t.ocm  tc  reac 
at  t .  oc«  s_per_  recc  ref 
,  nperreic 


i  oumoer  of  oytes  we  need  tc  put  ;rtc  the  o-'p-*  f:.e 
1-C  'the  size  cf  the  outp.t  vector  ;  f  we  hav*  -ere  * 


es  per  input  o.oc*  •  ..aoe.  <t 
process  •  tc* a  .  t . oc< s 
2  A  L ,  *ir :  n ; 


*nc  jg'  ir.p.t  data  t:  wr.te  c  .*  r..c:«*  a' 

'  '  •  ac,  se*  default  *  cut;.:  t.'  :«s  * ;  .  . *  terw.se 

‘t  e  tunoer  we  ca*  f.... 

:x  i  if  .oat  ii!  :r.  bl  2.  .  •  C . 1 
r  21,  oc. rc*s 

ad  •  : f : *  -  f.oat  tcta.  r.:c*s 
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«  float (blocka_to_read) )  ♦  0.5) 

3  -  1 
1  -  1 
k  -  1 

qult_flag  -  0 
qult”fiag2  -  C 
3rec  -  1 
raw  bytaa  -  0 

call  reader (blocks  to  read, bytes_per  block, disk, raw_bytes, 
&  qult_f lag, 3rec, Ibll, nblocks, qultj?lag2, input ^channel, 

t  diak_input_channel, blocks_per  set) 

call  transfar^status (aiabel, oldtlme, newtime, 
t  channels_to_process, number^  of_channals, process, 

t  hydrophones7fr»m*s_to_r*ad) 

if (quit_f lag2 .eq. 1 ) go”to  102 
if (quit_flag2 .aq.2) go  to  43 


C —  Now  throw  out  any  channels  we  don't  want  and  write  the  resulting  data 
C—  to  the  output  disk  file.  (watch  out  for  multiplexed  or  demuit iplexed 
C—  data; 


49 


43 


if (jrec.gt .frames  tojread)go  to  102 
1 f ( select or_mask  (£) . ne . C) then 

do  1-1, consacut iv*_channal_bytet 
output  data(j)  -  raw  byte_data{i) 

1  -  i  ♦  1 
3  -  J  ♦  1 

ifU.gt.raw  bytas)then 

if (quitJTlag.eq.l) go  to  43 
raw  bytes  -  C 

call  reader (blocks_tc_read, bytas  per_block, disk, 
raw_bytes, quit^f lag, }rec, iblk, nblocks, quit_flag2, 
input  channel, dlsk_input_channel, block s_per_set ) 

3rec  jrec  ♦  1  “ 

nrec  •  nrec  ♦  records  per^read 

if  (mod (nrec, 1C.) .aq. 0)then 

call  transfer_status (alabel, oldtime,  newtime, 
channel s_to_process, number_of_channels, process, 
hydrophones7frames_to_raadT  ~ 
end  if 
1  -  1 

if  (quit  flag2.eq.l)go  to  102 
if (quit^f lag2.eq. 2) go  to  43 
end  if 

if (}.gt .lfin) then 
oblocks  -  (3/512) 

call  output (disk_channel_out,  otype, oblocks, tblocks) 
if  (quit  flag.eq.ljgo  to  102 
3  -  1  ‘ 

end  if 
end  do 
else 

i  -  i  *  consacut ive_channel_bytes 
if  (i.gt . raw_bytes) than 
raw  bytes“-  0 

call  reader  (blocks_to_read,bytes_per_block, disk,  raw  bytes, 
quit_f lag, 3rec, iblk, nblocks, quit_flag2, input_channeT, 
dl BK_lnput_channel , block s_per_set ) 

3rec  -  3r*c  ♦  1 

nrec  -  nrec  ♦  records  per_read 
if  (mod (nrec, 10.) .aq. 7.0) then 

call  transfer_status (a label, oldtime, newtime, 
channel s_to_process, number_of_channels, process, 
hydrophone*7f rames_to_read) 
end  if 

if (quit_f lag2.aq. 1) go  to  102 
if (quit'f lag2.aq. 2) go  to  43 
i  •  1  -  raw_bytes 
if  (1 .gt . raw_bytes) go  to  50 
end  if 
end  if 

k  -  k  ♦  1 

if (k .gt .number_of_channels) k  -  l 


C- -  Update  program  status  on  terminal 
C 

go  to  49 


C —  If  wildcard  process  next  file  -  otherwise  let  user  choose  to 
C —  Process  another  or  exit.  Also  close  the  input  disk  file. 


Itime  -  ltime  ♦  1 
i f  (di *k .eq. 1 ) then 

call  dbclose  (disk_input_channel ) 
end  if 

i f  (ctype.eq. 'V* ) then 


C--  wos#  output  file,  but  first  calculate  the  number  of  records  in  the 
C--  f:.e  and  rewrite  the  header  block 

ilabel{3)  -  tbiocks  •  ilabel  (1)  /  ilabel  (2) 


a 
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call  dowrite  (disk  char.ne._out ,  raw_header,  1 ,  lerr,  O 
ca..  dbwait  (disK_cF.annei_out  i 
call  dDciose (diSK_channei_out ) 

e.se 

close { unit -2 ) 
enc  ;? 

if (wi idcard) then 
gc  tc  16 

else 

default  -  'K* 
another  -  getyr.  ( 

‘  Do  yo _  wart  to  process  another  file?’, 

Default  :s  to  exit  tne  program. .  • , 
default ) 

;f  ,  another)  go  to  19 
end  :  f 


tape  if  one  was  used. 


If  (dl sk  .  ne .  1 )  t her. 
write (5, •) ‘  Dismounting  Tape* 
mas*  *  amt  Sir  nounloaa 
status  -  sys^dlSTTvou  (cran,  %val  {mask;  i 
If:. not . status) call  .ids stop (%val (status) ) 
enc  If 
stop 


subroutine  check (iosb, iposi 


s  routine  checks  the  status  word  of  tne  lose  do  status  buffer) 
returnee  cy  system  ca.is  ana  aborts  the  program  If  problems  occur. 


lr.teger*2  1osd(4) 


If (lost (1) ,ne.  1)  then 

type*,'  error  lr.  system  cal*  at  posit  lor.  ‘ ,  Ipos 
do  1  -  1,4 

typelOC,  1,  iosb  (i ) 

format  (bx, • text_losb (•,  12,  ' )  -  ',18; 
enddo 
stop 
end!  f 
return 
end 


B2  -  Listing  of  Data  Output  Routine 


Subroutine  NAME:  OUTPUT 


Written  by: 
Joseph  B.  Farrell 
DREA 

15  Feb.  1986 


Latest  revision:  7  Aug.  1986 


subroutine  output (output_channel, filetype, oblocks, tblocus) 


This  routine  writes  sorted  output  data  into  a  disk  file  using 
the  DBMT  routines  from  NRL. 


output  data  (10240 


filetype 


integer* 4 
lnteger*4 
integer*4 


obiocks 

output_channei 

tblocks 


common  /odata/output_data 
tblocks  -  tblocks  ♦  oblocks 


C —  write  the  sorted  output  data  to  the  disk  file. 


if (f iletype.eq. *V' .or. filetype.eq. ’v* ) then 

call  dbwrite  (output_char.nel,  output  data,  obi  ocks ) 
call  dbwait (output_channel ) 

else 

wrlte(2) (output_data (lx) , lx-1 , (5l2*oblocks) ) 
end  if 
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return 

ena 


B3  ■  Listing  of  Mag-Tape  File  Open  Routine  "MTFELE" 


Subroutine  NAME:  MTFILE 


Written  by: 

D.  Vance  Crowe  i  Joseph  B.  Farrell 
DR  LA 

12  Mar.  1986 


Latest  revision:  15  Mar.  1986 


SUBROUTINE  MTFILE (input  channel,  header,  ifile,  wfile,  itime, 
verify, start_f fie, stop_fi le  ) 


Ihis  subroutine  opens  a  DREA  data  file  froiri  mao  tape. 


Farameter  and  internal  variable  declarations. 
IMPLICIT  INTEGER* <  (a-z) 

PARAMETER  MSG-5 


BYTE 

ICO 

! Null  byte  -  Null  character 

BYTE 

buffer (512) 

!Raw  header  data 

BYTE 

header (512) 

.'Raw  header  data 

CHARACTER *8C 

label 

IFile  name  labels. 

CHARACTER 

chan* 6 

!Used  in  nag  tape  channel 

CHARACTER 

default*l 

!User  input  default 

CHARACTER 

help (msQ)  *80 

!'Jser  help  prompts 

CHARACTER 

ifile*64 

! Input  filename 

CHARACTER 

lcnull 

! Null  byte  -  Null  character 

CHARACTER 

prompt  (msg)  *80 

!User  prompts 

CHARACTER 

tapejnark (10) 

! Checks  tape  marks 

CHARACTER 

tape  name (10) 

! Header  filename 

CHARACTER* 12 

VOLl“HDRl 

! VOL1 ,  HDR1  and  EOF1  labels 

CHARACTER 

wfile*  64 

! working  file  name 

INTEGER* 2 

f  ile_count 

INTEGER *2 

input  channel 

! Input  channel  for  tape 

INTEGER *2 

t«xt_Tosb (4) 

.’Status  buffer 

INTEGER-2 

lskip 

'skip  count 

INTEGER* 2 

start  file 

INTEGER-2 

stop_Tile 

INTEGER 

what_to_do 

'.option  flag  for  file  search 

LOGICAL 

end_of _tape 

! Ena  of  tape  indicator 

~OGICAL 

correct  file 

! File  flag 

LOGICAL 

getyr. 

iUser  response  function 

LOGICAL 

skip 

iSkip  fiie  flag 

LOGICAL 

verify 

equivalence  { 

tape  name (1 ), buffer (5) 

) 

equivalence  { 

label, buffer (1 )  } 

equivalence  ( 

icC, icnull  ) 

common  /io_st 

at  1  sties/ jrec, nrec 

external  ioS_skipf lie, lo$_readvblk, 

io$_skiprecord, io$_rewind 

Data  VCL1_KDR1 / ' V0L1KDR1E0F1 * /,  ICO/C/,  what _t o_do/-I / 
Data  for  user  interface. 


data  {prompt (i ), i-l,msg) / 

Enter  name  for  the  input  datafile.', 

Want  to  pass  over  this  fl le?  (-l-y/O-only/l-no/2-rewind) ' , 
Enter  a  name  for  the  Disk  output  file.'. 

How  many  data  frames  do  you  want  to  read?*, 

‘  Should  the  output  file  be  Vax  (V)  or  READRT  (R)  format?*/ 


data  (help (1 ), 1-1 ,msg) / 

4*  File  name:  UNIT  or  STRUCTURE : NAME . EXTENSION  * , 

4*  -1-pass  flieilook  at  next /C-correct /1-use  f 1 le/2-rewind . * 
4'  Default  name  will  be  the  same  as  the  Input  filename.*, 

4'  Frames  to  read  from  the  input  file.', 

4'  (V)  format  compatible  with  DISPVAX,  SASPEC,  etc.'/ 


data  1  file/ ' msaC : q38334 . dat '/ 
chan-'  '//ifile  (1:5) 


APPENDIX  B 


N-*. 


•.vj 


I 


>  .v 


e 


.VV> 

\<s 


APPENDIX  B 


i  i  [  i  t  ime .  eq .  I )  t  he  r. 
fiie_count  -  C 
call  gettap(chan) 


Assign  a  channel  to  the  mag  tape 

stat  us»«sysSassign  (char.,  input_channei , , ) 
if  i. not .status) call  1 ib$sccpT%val (status) ) 
ena  if 

L  i  (fi  ie  count .  gt .  stcp_file)  stop  'All  flies  processed.1 


Read  the  R7-11  heaaer  (incl.  filename;  from  tape  into  'buffer* 
i  error  -  0 

status-sysSqiow  (,  %val  (ir.put^cnannei } ,  ioS_readvbik, text_iosb, , 

, %ref (buffer) , %vai (512) , ,7# > 
if  i. not. status) call  libSstcp (%val (status) ) 

if  { .not .status. or. text_iosb (2) .eq.C) then  : t here  was  an  error 
ierror  -  1 

if  (text  lose (2) . ne . C) then  : there  was  a  tape  read  error 
else  ! cr  there  was  ar.  EOF 

end  if 
end  if 

If ( ierror. Eq.l)  Go  to  11C  ! Try  again 

itype  -  INDEX {  VOLl^HDRl,  label (1:3)) 

1 f { i type . Eq . 0 ) then  " 

If  {  itime.eq.i  jthen 

1  skip  -  -2  (Probably  a  data  record  -  go  bac*  to  file  st 

e.se 

isxip  “  ♦ 2  .’Probably  a  data  record  -  gc  ahead  tc  file  s 

end  if 
end  if 

I f ( itype .  Eq .  1  )  Go  to  110  1  VOL  label,  read  next  header 

I f ( itype . Eq . 5  )  Go  to  150  ! KDR  label,  checx  names  etc 

I f (itype. Eg. 9  )  iskip  -  +1  'EOF  label,  skip  forward  1  EOF 

!  Shipping  file  rr.arics 

status-sysSqiow  (,  %vai  (ir.put_channei ) ,  io$_skipfiie, 

1  text_iosb, ,, %val (iskip; ) 

if  ( .not . status) call  iibSstop (%val (status) ) 
call  check  (text_iosb, -3) 

Go  to  11C 


Check  tc  see  if  we've  reached  the  end  of  this  tape  (Filename  all  0's) 

end  of  tape  -  .true, 
do  I  -“1 , 1C 

if  (tape  name(i)  .ne.  ienuil)  then 
er.c_o?_tape  -  .false, 
endi  f 
end  do 

if  (end_of_tape)  then 

type*7'End  of  tape  encountered,  Rewinding.* 

Rewind  tape ! . . . 

status-sysSqiow  {, %val  { input _c nan nt. ) , ioS_rewind, text_icso,  ,,,,,,,) 
i f ( . not .status) call  1 ibSstop ■ %va.  (stat us) ) 
wnat_tc  do  -  -1 

Gc  to  lie  !Penu  the  volume  .a be Is 

end  if 


If  not  at  end  of  tne  tape,  cnec*  *.c  see  if  we're  picking  a 
numbered  file.  I f  we  are,  pos.titr  tne  tape  tc  that  file. 

if  (start_file  .gt .  1  .and.  it  ime.eq.  I .  t :  er. 
iskip  -  3  •  (start_file  -  1; 
file_count  -  1 

type* , 1  Skipping  to  next  f:.e.’ 

Go  to  1 20 
end  if 

if  (start  f  i .  e .  eq .  1  .and .  it  i~c- .  eq  Is:.  ie_  count  - 


See  if  we're  reading  tne  c crre~-  «  .  (wilccaras  allowed; 

correct  file  •  .true. 


If  we're  doing  numbered  file",  w«  cc-'t  .have  tc  check  the  filename. 

if  (start_f  i  le.r.e  .  C  gc  it 
*  -  index  (i file,  ' : ’ :  ■  1 
period_f lag  -  C 
name  length  •  .er.2  :f..‘- 


C.neCK  for  wildcard  '**  c:.a:<>  nput  fi.ename. 


if  ( i  2 .  re .  A : ;  t  r.e 

Compare  input  ar.d  tape  fiie\n"*  :  :ter  if  .  np.,1  ret  1  •'. 
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11-lchar  (tape  named)) 

If  dl.ne.12)  then 

correct  file  -  .false, 
go  to  21 
end  If 
else 


If  we  got  a  wildcard  character,  look  for  ■ .•  If  period_flag 

If (period  flag.eq.O)then 

do  while  (lchar (tape_name (1) ) .ne. 46. and. 
i.lt.10) 

1-1  +  1 
end  do 
1-1-1 
peri od_f lag  »  1 
else  - 


If  we've  already  passed  the 


then  filename  is  assumed  correct. 


go  to  20 
end  If 
end  if 
3-1  +  1 

if (3.gt.name_length)go  to  20 
end  do 

If  this  was  the  correct  file,  indicate  success  and  proceed. 

If  (  correct  file)  then 
do  k  -  1, Index (ifile, •:' ) 

wflle(k:k)  -  lflle(kik) 
end  do 
do  k  -  1,10 

j-k+  index (iflle, ':•) 
wflle(j:j)  -  tape_name (k) 
end  do 

type*,'  file  1 , wflle (6: 15) , 1  found' 
goto  220 

end  if 

If  it's  not  the  right  file,  skip  to  the  next  one  (giving  the  user 
the  option  of  processing  the  file  which  has  been  found) . 

continue 
if (verlfy)then 

type*,'  Found  file:  ',tape_name 
If (what  to  do.Eq.  -1)  Then  !Ask  the  question  the  first  time 
call  input! (prompt (2) , what  to  do,help(2)) 

Endlf  ~  ~ 


what  to  do  -  0 


If (what_to_do  .Eq.  0)  Go  to  180  IFind  the  correct  file 

If (what~to~do  .Eq.  1)  Go  to  190  !0se  this  file  as  correct  one 

If  (what_to_do  .Eq.  2)  Go  to  160  JRewind  tape. 

If <what_to_do  .Ne.-l)  Then 

Type  *,'  Assume  you  want  to  skip  to  next  file 
endlf 
iskip  -  +3 

type*,'  Skipping  to  next  file.' 

Go  to  120  .'Skip  forward  3  f  1 1 


.'Skip  forward  3  file  marks. 


The  user  wants  this  file  that  he/she  found. 

do  k  -  1, index (iflle, ) 
wflle(k:k)  -  ifile(k:k) 
end  do 
do  k  -  1,10 

3-k+  index (iflle, ':') 
wflle(j:J)  -  tape_name (k) 
end  do 
go  to  220 


Skip  over  file  mark  to  the  beginning  of  the  drea  header  data 

status-sysSqiow (, %val (lnput_channel) , io$_skiprecord, 

1  text  iosb, ,  ,»val  (1) , , ,, , )~ 

if (.notlstatuslcall  llbSstop (aval (status) ) 

We're  finally  in  position,  so  read  in  the  DBEA  header 

status-sysSqiow (, aval ( input _channel) , io$_readvblk, text_iosb, , 
1  ,tref (buffer), %val (512),,,,! 

if (.not.status)call  1 ibSstop (%val (status) ) 

Type  *, '  finished  reading  drea  header' 
call  check  (text  iosb, -5) 


At  this  point,  we've  read  in  the  header  from  tape. 

file  count  -  flle_count  ♦  1 
DO  3!C  11-1,512 
header (1 1) -buffer (11 ) 
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RETURN 

END 


B4  -  Routine  to  Allocate  and  Mount  Mag-Tape 


c** . . . . . . 

C  Subroutine  NAME:  GETTAP 

C**- . . . . 


c 

c 

c 

c 


Written  by: 
Joseph  B.  Farrell 
DREA 

25  M jr.  1986 


Latest  revision:  7  Aug.  1986 


subroutine  get tap (chan) 


structure  /ltmlst/ 
union 
map 

INTEGER* 2  buflen 

INTEGER-2  code 

INTEGER-4  bufadr 

INTEGER *4  end 1st 


end  map 
map 

INTEGER-4  end_list 

end  map 
end  union 
end  structure 


record  /ltmlst/  mnt_llst  (3) 
include  *  (Smntdef) 

Include  ‘ (Sdmtdef) ' 

Include  •  (Sssdef } * 


CHARACTER 
INTEGER -4 
INTEGER-4 
INTEGER- 4 
INTEGER -4 


chan -6 
mask  /C/ 
status 

sysSmount 

sysSailoc 


Used  ir.  mag  tape  u  ha  fine 


C 

C 

C- 


Allocate  and  mount  the  mag  tape, 
status  -  sysSailoc  (char., , , , ) 


If  tape  is  already  mounted  snip  mounting  section. 

if  (st at us. eg. ss$_devmount ) then 

type*,’Tape  is  already  mounted  -  assuming  yo-  did  it.1 
go  to  1C 
end  If 

if  (status.eq.  ssS_deval  *ocj  tner. 

write  (5,-)*  Error  -  device  pro&ac.y  a.,ccat»c  tc  another  user' 
stop 
end  if 

if  (status. eg . ssS  devalralloc: ther 

write (5, • :  *  Device  already  allocated  to  yo«. * 
end  if 

mask  -  C 

mask  •  mnt$m_foreigr.  .or.  mnt  Sir  _mes  sage  .or.  mr.t  Sr  noaian*. 

mnt_list (1 > .buflen  -  6 
mnt~list  (1)  .  code  -  mnt S  devnar 
mnt_iist  (1  j  .bufadr  -  Hoc  (chan: 
mnt~li«t  (1 )  .endlst  -  C 
mnt_llst (2) .buflen  -  4 
mnt_list  (2) .code  •  mntS^flags 
mnt~list (2) .bufadr  -  %loc(mask> 
mnt”list (2) .endlst  -  C 
mnt“l 1st (3) .end_i 1st  -  C 

status  •  sy sSmour.t  (wni_l lit 

i f ( . not . stat us) cal .  1 ibSstop ‘%va .  states  ■ 

continue 

return 

end 


B5  -  Channel  Selection  Subroutine 
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Written  by: 

Joseph  B.  Farrell 
DR£A 

21  Feb.  1906 

Latest  revision:  6  Jui.  1906 

This  subroutine  chooses  which  of  the  channels  ir.  an  RT-11  or  VAX 
format  oata  file  are  to  oe  processed.  It  uses  the  GETVEC  subroutine 
written  by  Doug  Peters  (ASP  Sumer  Student)  to  read  a  vector  iron,  the 
user’s  tenrlnal. 

SUBROUTINE  CALLS: 

1)  INPUTI  -  Routine  which  reads  ar.  integer  from  the 


)  INPUTI  -  Routine  which  reads  ar.  integer  frotr  the 

terminal 

2)  GETVEC  -  Routine  which  reads  a  vector  frotr  the 

texir.lna* . 

3:  WAIT  -  Waits  for  a  specified  amount  of  time. 

4:  £rase_screer.  -  Clears  the  terminal,  screen. 

t.  Set _curscr  -  Moves  the  cursor  tc  a  specified  location 


subroutine  channel _se*ect (blabei, number  of^cnanneis, out spec, itime, 
4  channe.s_ tc  process,  process, nySrophonea, si ft_f lag, wfl le, 

4  wi.ocaro7sKlp) 


Parameter  ir.a  internal  variable  declarations. 


implicit  integer* 4  (a-*) 


parameter  msg-1 


character 

character 

character 

character 

character 

character 

integer *2 
integer*2 
integer *2 
integer"! 
integer*2 
int eger • 2 
integer*2 
integer* 2 
integer*2 
Integer *2 

integer*4 
1  nt eger  *4 
integer* 4 
integer *4 
integer *4 
int eger *4 

LOGICAL 

LOGICAL 

LOGICAL 

LOGICAL 

LOGICAL 

LOGICAL 


blaoel (120) 

default *1 
help  irosg) •? C 
wfiie*64 
prompt (mag) *?C 
response* 3 
string-12 

channels (120) 
orlgina*_channels  (120) 
r:ydropnon0»  (22$) 
new  HP 

nuaEer_of _channe 1 s 
number  of  hydrophones 
phones7l2ff) 
original  phones(120) 
process (T20) 
sift_flag 

channel s_to_pro cess 

origlnal~channels_to_process 

cnees  channel  ~  * 

f  ir.isT) 

if.ag 

ip 

all  flag 
ok (120) 
out  spec 

preserve_phones 

proceed 

wildcard 


Lata  for  user  interface, 
data  (prompt ( i i ,!•! ,msg) / 

4'  How  many  channels  do  you  want  to  process?  (-1  for  all)*/ 

data  (help (1) , i-l,KSC) / 

41  Enter  the  number  of  channels  to  process.'/ 

Variable  l r.it  lal  i sat  iona 


number_of  hydrophones  -  C 
new  HP  -  Z 

i f (Tt ime.eq . 1 .or . outspec) al l_f lag  -  .FALSE. 


Letermine  which  channels  contain  acoustic  data.  Collect  the 
rydrophone  numbers  sna  corresponding  channels. 


dc  !• 1,2 ‘number  cf  channels, 2 
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K  -  j/2 

if  (blabel (i)  .ge.  Other 

number  of_hydropnones  -  n umber _of_hydrophones  ♦  1 
channel* (number  of_hydropnones)  -  k 
phones (number_o7_hydrophones>  -  blabel (i) 
end  if 
end  do 


Must  be  at  least  one  hydrophone  in  the  file. 

if  (number_of  hydrophones  .It.  1) 

4  stop  ' ChanneT_select  -  no  acoustic  data  found  in  file.' 


Allow  channel  numbers  or  H/P  numbers  to  be  preserved  between 
files. 

1  f (.not.outspec.and.itime.eq.l . and. wildcard) then 
default  -  ‘Y* 
preserve_phones  •  getyn( 

4  '  Select  the  same  H/P  (Y)  or  channels  (N)  from  each  file?', 

4  *  Default  will  select  tne  same  H/P  from  each  input  file.', 

4  default) 

end  if 


If  itime  >  1  compare  hydrophones  with  those  from  the  original  file 
and  write  out  a  warning  if  the  chosen  ones  differ. 

if ( .not .outspec.and.it ime.gt . 1) then 
misting  f lag  -  C 
1  -  1 


Loos  tnrougn  the  available  pnones  to  see  if  the  desired  ones  are 
there . 

do  j  ■  1,  orlginal__channels_to_process 
do  k  -  l#number”of_hydrophones 

If  (phones  (s) . eq. origins l_phones ()) ) then 


Use  the  process  vector  to  point  at  the  location  of  the  desired 
H/P  if  we're  keying  on  H/P. 

if (preserve  phones) then 
process (I)  -  k 
1-1*1 
end  if 
go  to  1 
end  if 
end  do 


Set  a  flag  to  Indicate  H/P  missing  if  we  can't  find  one. 

missing^flag  -  1 
continue 
end  do 


If  there  are  H/P  missing  write  a  warning  and  let  the  user  decide 
wnat  to  do. 


if (missing^flag.eq. 1) then 
call  erase_screen (1 , 1 ) 
call  set_cursor (2, 1) 

call  vtmess ( ' re ' , 'c* ,  '  Channel  Setup  ’) 
type  1000,wfile 

1CCC  formate  Working  on  file:  '  ,A20) 

type*,'  The  h/p  available  differ  from  the  originals.' 
type  *# •  OLD:  • , (original_phones (} j) , j j»l, 

4  original_channels_to  process) 

type  *7*  NEW;  (phones (jj) 

4  number_of_hydrophones) 

if ( . not .preserve_phones) then 

type*,'  You  are  keying  on  channels  rather  than  phones  ' , 
4  'so  I  am  proceeding.' 

call  wait  CO  :  :3  ' ,  5) 
end  if 


C- 

C 


If  we're  preserving  a  set  of  H/P  betwen  files, 

Option  to  proceed  with  the  found  subset  of  H/P  or  to  sxip  file. 

If (preserve_phones) then 
proceed  -  .false, 
if  (1 .gt .1) then 
default  -  'Y* 
proceed  -  getyn ( 

4  'Proceed  using  the  subset  of  requested  phones  found?', 

4  'Re-speclfy  H/P  or  skip  this  file  if  reply  :s  K.'r 

&  default) 

end  if 


If  the  user  chose  not  to  proceed  with  the  subset  found,  :•  ve 
the  option  on  re-specifying  the  H/P  or  of  skipping  thf  f  e. 
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1 f ( . not . proceed ) t  hen 
default  -  *Y* 
proceed  •  getyn( 

'Skip  to  the  next  flic  in  the  input  set?', 

‘Default  is  to  request  •  now  sot  o f  H/P  to  proooas.’, 
default) 
if (proceed) then 
akip  *  1 
rot  urn 

ol  so 

now_HP  -  1 
go  to  8 
ond  if 
ond  if 
ond  if 


If  wo* re  using  tno  subsot,  continuo  processing, 
ond  if 


If  wo'ro  keying  on  cnanneis  rather  than  H/P  make  sure  wo  have 
enough  channels. 

If (. not .preserve_phones) thon 

Only  do  it  if  wo'ro  not  processing  all  channels  in  the  file. 


if (.not .all_f lag) thon 
k  -  1 

do  j  -  1, origins 2 _channels_to_proces# 


Cneck  to  so*  if  a  requested  channel  is  >  than  the  #  wo  have. 

if (orlginal_channels (}) ,gt .number_of_hydrophones) thon 
type*,*“  You  asked  for  a  channel~not  found  in  the*, 
’  input  file. ' 

If  channel  out  of  range,  user  can  skip  file  or  proceed. 


default  -  ' Y ' 
proceed  -  getyn  ( 

4  ‘Snip  to  the  next  file  in  the  input  set?*, 

4  ‘Otherwise  use  this  file  with  reduced  t  of  chans.', 

4  default) 

if (proceed) then 
akip  -  1 
return 
end  if 
else 

type*,'  Proceding  with  reduced  •  of  channels.' 
process(k)  -  original  channels(j) 
k  -  k  ♦  1 
end  if 
end  do 

channels_to  process  -  k  -  1 
end  if 
end  if 
end  if 

i f (all_f lag) then 

channel s_to ^process  -  number_of_channels 
go  to  5 
end  if 

c . 

C —  Clear  the  screen  and  write  a  header  message, 

c 

if (ltime.eq.l.or.outspec)then 
8  call  erase_screen (1, 1) 

call  set_cursor (2, 1) 

call  vtmess ('re*, 'c', •  Channel  Setup  ') 

c . 

C -  Show  the  acoustic  channels  available  and  let  the  user  choose  some. 


type  1,  number  of  channels 

formate  ',t3,T  Tfvere  are', 13,'  channels  in  the  input  file.’) 


2 


4 


3 


4 


type  2,  (channels (i) , i-1, number_of_hydrophones) 
formate  ',t3,'  The  following  are  acoustic  channels:  '  / 

■  ',t3, 26(12:  \  ')/'  ',t3, 26(12:  \  ')/*  ■ , t3, 26  (12: • , • ) ) 

type  3,  (phones (1) , i-1, number  of_hydrophones) 
formate  ',t3,'  The  correspon3lng  H/P  numbers  are:  '/ 

'  ', t3, 26(12: ',') /'  ', t3, 26(12: ', ')/'  ' , t3, 26 ( 12 :  • , * ) ) 


channel s_to  process  -  -1 

A  call  input lTprompt (1) , channels  to_process, help (1 ) ) 

end  if 

C . 

C-—  Set  flag  Indicating  all  channels  to  be  processed  for  subsequent 
C—  files  if  necessary. 

if (It lme.eq.l .or . new_HP .eq. 1 .or. out spec) then 
1 f ( channel s_t  oprocess . eq . - 1 ) then 
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ailjflag  -  .TRUE. 

cnannela  to_procesa  -  nuaber_of_channel * 
end  If  ~ 

orlginsi_channela_to_process  -  channels_to_process 
•  no  if 

User  must  choose  hydrophones  >-  1  and  <-  number  available 

if  (channels_to_process  .It.  1  .or, 

channel s_to_process  .gt.  number_of_channela)  than 
type  3500 

format!/*  channei_select  -  illegal  number  of  channel*  to*, 

’  process  chosen.'/) 
call  wait ( *C  ::2*,5) 
call  erase_screen (5, 1 ) 
go  to  4 

Load  all  H/?#  channels  and  gains  if  all  are  to  be  analyzed. 

else  if  ( channel s_to_process  .eg.  numbered ^channels)  then 
do  4CC  i-1, number_of_channeis 
process(i)  -  i  ” 
hydrophones {1 )  -  phones  (i> 
if  (ltime.eq. Dthen 

original_phones  (i)  *  phones (i) 
originai_channeis (i )  -  process (i) 
end  i  f 
continue 
return 
end  if 

Identify  individual  h/P  if  a  subset  of  the  total  was  chosen. 

if (itime.eq.l.or.new  HP .eq. 1 .or .out spec) then 
type  400C 

format (/'  ',t3, '  Enter  the  cnannels  you  wish  to  study  :  *,$) 
call  getvec (val, ok, channel *_to_process) 
if  (ok  (1) )  then 

do  )-l, channel s_to_process 
process (3) -int  (val  (3) ) 
end  do 

Set  flag  to  indicate  that  a  channel (s)  must  be  removed  from  data. 

si f t_f lag  -  1 

else 

type  '(/••  channel  select  -  sorry,  no  default. ••)* 
call  wait  CO  :  :2* ,  5) 
call  erase  screen (5,2) 
go  to  4 
end  if 
end  if 

do  3-1, channels  to_process 

hydrophones (3)  *  phones (process (j) ) 
if (itime.eq. 1) then 

original_phones (j)  -  phones (process (3) ) 
original_channels (3)  -  process (3) 
end  if 
end  do 


B6  -  Routine  For  Choosing  Data  Segment  To  Be  Transferred 


SUBROUTINE  NAME:  SKIPPER 


Written  by: 

Joseph  B.  Farrell 
DREA 

14  Feb.  1986 

Latest  revision:  23  May.  1986 


This  program  uses  the  header  to  determine  the  starting  time  cf  ar, 
input  file,  asks  the  user  to  enter  a  desired  start  time,  and  then  steps 
into  the  file  the  desired  amount. 


subroutine  skipper (lnput_char, points_per_block, samp! lng  frequency, 
1  alabel, blocks_per_record, oldtime, newtime, disk,  diak_start  block, 

1  file_typ#,outspec7itime, total_biocks, block_t ime, frame_tlme, 

1  frset ional_biocks, number_of_channels) 

Parameter  end  internal  variable  declarations. 
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IMPLICIT  INTEGER-4  (a-l) 


PARAMETER  asg-2 

CHARACTER  alabal *32C 

CHARACTER  default *1 

CHARACTER  heip  (msg)  •  "7C 

CHARACTER  newtime-8 

CHARACTER  oldtime-8 

CHARACTER  transtl  nm-8 

CHARACTER  prompt (msg) *70 

INTEGER  blocks_per  record 

INTEGER  blocks  to  Bklp 

INTEGER-4  disk 

INTEGER-4  disk  stsrt_block 

INTEGER-2  flle~type 

INTEGER-2  input_chan 

INTEGER-4  itime 

INTEGER-2  new  tlme(3) 

INTEGER *2  number_of_channeis 

INTEGER -2  redo 

INTEGER-2  text  losb{4) 

INTEGER* 4  total  blocks 

INTEGER-2  trans?ex_t ime (3) 

LOGICAL  getyn 

LOGICAL  newtim 

LOGICAL  outspsc 

LOGICAL  Status 

REAL  begln_tlme 

REAL  block  time 

REAL  f  rsctlonal_blocjcs 

REAL  frame  time 

REAL  old_tTroa {3) 

REAL  points__per_block 

REAL  block a_per~hour 

REAL  raal  skip 

REAL  sampling  frequency 

REAL  tapa  tins 

REAL  tast 

REAL  t irae_difference 

EXTERNAL  io5_skipf 1  la, io$_raadvblk, io$_skiprecord, io$_rewlnd 

Data  for  usar  lntarfaca. 

data  {prompt (1) , 1*1, msg) / 

4 '  Do  you  want  to  bagin  procassing  at  soma  othar  tin*.*, 

4*  Enter  tha  tlma  at  which  you  wish  to  start.*/ 

data  (halp (1) , 1*1, msg) / 

4'  Dafault  Is  to  start  at  tha  tlm*  shown.*, 

4*  Format  Is  HH:**i:SS.'/ 

Variable  Initializations. 

rawtiin  -  .falsa. 
oldtimaU:)  -  alabal  (12:19) 
nawtimad:)  -  alabal  (12: 19) 
if  (file_type.eq.l)then 

block s_per_hour  •  3600. *sampllng_fraquancy/ 

1  (point s_par_block) 

block  tlma  -  1. /blocks  par  hour 
and  If 

Print  tha  file  start  tlma  on  tha  user*s  terminal. 

newtim  -  .FALSE, 
if (ltime.eq. 1 .or .out spec) than 
call  erase_screen (1, 1) 
call  set_cursor (2, 1) 

call  vtmess { 're* , 'c' , '  This  file  starts  at  * //alabal (12: 19) , • 

Let  tha  usar  determine  a  start  tlma  for  data  analysis. 

(If  NEWTIM  Is  returned  as  "FALSE"  analysis  starts  at  tha  beginning 
of  tha  file.) 

dafault  -  'N' 

newtim  -  gatyn (prompt (1) , halp (1) , dafault ) 
and  if 

lf(  newtim  )then 

Decode  tha  file  start  tima  into  tha  vector  Oid_time. 

decode (2, 1000, alabal (12: 13) ) old  tima (1) 
decode (2, 1 00C, alabal (15: 16) > old‘t ime (2) 
decode (2, 10CG, alabal (16:19) )old~time (3) 
format (f 2.0) 
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:  f  we're  chang; ng  start  times,  reac  ;r.  tne  new  tin*  ana  "DECODE-  ; 
;r.tc  the  array  New_tlme. 

cal.  1 nput  a (prompt { 2 ) , newt ime , ha 1 p ( 2 ) ) 
i f (newt ime. eq.oidtime) go  tc  1 
cecode (2, 1001, newtime (1:2) ) new_t ime (1) 
aacooa  (2, 1001, newt iroe (4:5)) new_t ime (2) 
dacoda (2, 10C1, newt ima (7:8) ) new~t ime (3) 
format  (12) 


Tape_tlme  is  mat  to  tna  input  file  atari  tlroa  ilr.  aacimal  nours) . 
&eg;r_tima  ia  aat  tc  the  procaaaing  start  time. 

tape  tima  -  old_tima(l)  ♦  (old_tlme(2)  /  60. C)  ♦ 

1  (oi3_time(3)  /  360C.0) 

oagin^tlma  -  f loat (new^t ima (1 ) ) ♦ (f ioat <new_t ima (2) ) /  60.0)4 
(float (new_t ime (3) )  /“360C.C) 
t ime_dif ferenca  -  begin_time  -  t ape_t lme 


Ka*e  sure  we  start  on  a  block  which  oegine  with  the  first  channel 
1!  we're  acing  a  .DAT  file. 

nDlocks  -  1 

if  (fract lcnaljolocks.gt .0.0001) than 
do  nPlocaa  -  2, number  cf__channels 

teat  •  float  (noloclcs^fract lonal  blocks 

if{  laba(teat)  -  abs  (int  (teat) ) )  7lt  .O.OOODther. 

go  to  22 
and  if 
and  do 
anc  if 


Cnee*  tc  maice  aura  we  aren't  moving  into  tna  middle  of  an  FFT  or 
a  spectrum  if  the  file  ia  . FTR  or  .  PWR 

if (file_type.eq.2) call  timer (t lme^dif ference, frame_time, redo) 
if (fiie~type.eq.3) call  t imer (t ime~dif ference, frame  time, redo) 
if (redc.eq.l) go  to  21 


Determine  how  many  records  to  skip  before  beginning  processing. 
(fiie_type-l  indicates  .DAT,  -2  indicates  .FTR,  and  -3  .PWR) 

if  (file  type.eq. 1) then 

biocEa_tc_aklp  -  int (blocks_per_hour  •  t lme^dif ference) 


Make  sure  we  start  on  a  block  which  begins  with  the  first  channel 
if  we're  doing  a  .DAT  file. 

if (amod (float (blocks_to_skip) , float (nblocks) ) .eq. 
t  0.)ther 

go  to  23 

else 

blocks  to_skip  -  blocks_to_skip  -  1 
go  to  74  ~~ 
end  if 

real_sklp-tape_tlme  ♦  (float  (blocks  to  skip) /blocks  per  hour) 
else  ~  ~ 

blocka_to_sklp  -  (int (t lme_dlf ference  /  block  time)) 
real_8Kip-tape_time  ♦  (f *oat (blocks_to_skip)  •  block  time) 
ena  if 


Set  the  start  time  to  the  time  we're  actually  going  to  skip 
*nto  the  file  (We  may  not  be  able  to  skip  exactly  to  the  requested 
time  because  of  the  finite  record  lenght  in  the  input  file). 

new  timed)  *  int  (real_skip) 
real_skip  -  (reol_skip-new_t lme (1) ) *6C 
new  time  (2)  -  intlreal_skip) 
reaT_skip  -  (real_skip-new_t ime (2) ) *60 
new_tlme{3)  -  nint (reai_skip) 


"ENCODE"  the  actual  start  time  into  the  character  stinQ  NEWT IKE 
and  display  It  or  the  user's  terminal. 

encode (2, 1001 , newtime (1:2))  new_time (1) 
encode (2, 10CZ, newtime (4:5) )  new_time (2) 
encode (2, 1001, newt lme (7:8) )  new  time (3) 
type*, '  Actual  start  time  wili'be  newtime 
a label (12:19) -newtime (1:8) 


Stip  the  requested  number  of  records 
. . (dl ak.ne. 1 ) then 


S*.ip  recoroa  or  tape  if  that  ia  the  media  being  used. 

If  (blocks_to_akip  ,ne.  0)  then 

type*, blocks  to_skip, '  physical  blocks  will  be  skipped.' 
atatua-ayaJqlow {, Ival  (input  chan>,io$  akiprecord, text  lost, 
1  , 6va a  (biocns_to_akip) ,,,,,) ~ 

If  (.not .status) call  llbSatop  (%val (status) ) 
cal.  check <text_iosb, -11 ) 
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•ndif 

tilt 

c . 

C -  Sat  the  start  block  for  reads  from  disk. 

C 

disk  start  block  -  blocks_to  skip  +  1 
type7,'  Disk  start  block~wlTl  be:  • , disk_start_block 
•nd  If 

call  wait  CO  ::2\5) 

•  ls« 

C . 

C —  Control  jumps  here  if  we're  going  to  start  at  the  baginning  of  the 

C  fils, 

c 

1  disk_start_Mock  -  1 

•nd  if  ~ 


C - 

C 


20 


C.  . 

c— 

c 


Now  give  th«  user  th«  option  of  specifying  the  number  of  records 
to  be  read  from  the  input  file  or  specif lying  a  time  Interval. 

if (itime.eq.l.or .outspec) then 
total  blocks  -  -1 
call  Input i ( 

4  *  Enter  n  to  X-fer  n  blocks, -1  for  all, -2  to  specify  tints*, 

4  total_blocks, 

4  *  -2  will  let  you  enter  a  time  interval  for  the  transfer.') 

if (total_blocks.eq.-l) then 
total  blocks  -  100000 

else  “ 

if  (total_blocks.eq.-2) then 
transtime  -  '00:01:00* 
call  Inputs ( 

4  '  Enter  the  length  of  time  of  the  transfer.’ , 

4  transtime, 

4  '  Format  is  HHsMCSS.1) 

decode (2, 1001, transtime (1:2)) transfer_t ime (1) 
decode (2, 1001, transtime (4:5) ) transfer_time (2) 
decode (2, 1001, transtime (7:8) ) trans foretime (3) 
time_dlfference  •  transfer  time(l)  ♦  Itransfer  time  (2)  /  60.0)  ♦ 
1  (transfer  time (3)  /  3600.5) 

if (time_Hif ference.eq.0) then 

type*,'  You  specified  a  zero-length  transfer  -  try  again.' 
go  to  20 
end  if 


Check  to  make  sure  we’re  transferring  at  least  a  full  FFT  or  spectrum 

if (file_type.eq.2) call  timer (tlme_dif ference, frame_time, redo) 
if (file_type.eg.3) call  timer  (time_dif ference, frame_time, redo) 
if (redoTeq.l)go  to  20 

total_blocks  •  time  difference  /  block_time 
if (total_blocks.eq.C) then 

type*7‘  WARNING  -  you  are  trying  to  transfer  0  blocks!', 
4  '  Try  again. ' 

type*,'  Transfer  time  must  be  at  least', 

4  frame  time*3600. 

go  to  'JO 
•nd  if 
else 

total_blocks  -  total_blocks  •  block s_per_record 
end  if 
end  if 
end  if 
ret  urn 
•nd 


subroutine  timer (time_dif ference, unit_time,  redo) 

IMPLICIT  INTEGER* 4  (a-r) 

INTEGER* 2  choose 

INTEGER*2  redo 

REAL 
REAL 
REAL 
REAL 

redo  •  0 

check  -  amod  (tlme_di f ference, uni t_t ime) 
seconds  -  unit_tlme  *  3600. 
if ( check. ne.0. It hen 

type*,'  You  are  not  using  an  integral  number  of  records.' 
choose  -  1 
call  input i( 

4  '  Type  1  to  take  closest  record  start,  2  to  re-speclfy  time', 

4  choose, 

4  '  Closest  may  be  earlier  or  later  than  the  chosen  time.') 

if (choose. eq. 1) then 

time_dlf ference  -  anint (time_dif ference  /  unit_time)  • 

4  unit  time 

else 


check 

seconds 

t ime_d 1 f f er • n ce 
unit  time 
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type*,  1  Tim*  difference  must  De  a  multi.pl*  of, seconds, 
'  seconds . ' 

•  redo  -  1 
•nd  If 
•nd  if 
return 
end 


B7  -  Routine  to  Decode  Switches  in  Filename 


SUBROUTINE  NAME:  SWITCHES 


Written  b> : 
Joseph  B.  Farrell 
DREA 

3  Jui.  1986 

Latest  revision:  3  Jui. 

1986 

This  subroutine  picks  switches  from 

a  user  input  filename. 

subroutine 

switches  (  filename, verify, start_flle, stop_file) 

character 

filename-64 

INTEGER-2 

number 

INTEGER-2 

start  file 

INTEGER-2 

stop  file 

INTEGER-2 

start 

INTEGER-2 

stop 

logical 

verify 

Initialize  variables 


heck  to  see  if  verify  flag  is  present 

start  -  index  (  filename  ,  '/V') 
lf{  start  .ne.  0  )  verify  -  .TRUE. 


Look  for  a  START  flag  ,  and  if  present  decode  the  starting  file 
number. 

start  -  index  (  filename  ,  ‘/START-*  ) 
if  {  start  ,ne.  0  )  then 
start  -  start  ♦  7 

stop  -  index  (  filename (start : )  ,  •/*  ) 
if  (  stop  .eq.  0  )  then 
stop  -  len2  (filename) 
number  -  stop  -  (start-1) 

else 

stop  -  start  +  stop  -  2 
number  -  stop  -  (start-1) 
end  if 

decode  (number, 100, filename (start : stop) )  start  file 
format (  12  ) 


Look  for  a  STOP  flag  ,  and  if  present  decode  the  stopping  file 
number. 

start  -  index  (  filename  ,  '/STOP-*  ) 
if  (  start  .ne.  0  )  then 
start  -  start  ♦  6 

stop  -  index  (  filename (start : )  ,  '  *  ) 
if  (  stop  .eq.  0  )  then 
stop  -  len2 (filename) 
number  -  stop  -  (start  -  1) 

else 

stop  -  start  ♦  stop  -  2 
number  -  stop  -  (start  -  1) 
end  if 

decode  (number, 100, filename (start : stop) )  stop  file 
else  ~ 

stop  file  -  start  file 
end  if 
end  if 


Look  for  switch  mark  and  remove  all  switches  from  the  filename. 


start  •  index  (  filename  ,  V) 
if  (  start  .ne.  0  )  then 

filename (start :len2 (filename) )  • 
end  if 
return 
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end 


B8  •  Routine  to  Determine  Wildcard  File  List 


Subroutine  NAME:  Disk  Wildcard 


C*  "  * 

c 

c 

Written  by: 

c 

Josaph  B.  Farrail 

c 

DREA 

c 

10  Jul.  1986 

C 

Latest  revision:  1  Aug.  1986 

subroutine 

dlsjc^wildcard  (iflle,  ltlme,  wfile) 

CHARACTER 

string* 80 

CHARACTER 

1  file*  64 

CHARACTER 

wfile*  64 

INTEGER-4 

dev  flag 

r 

INTEGER-4 

itime 

Do  a  directory  using  the  input  filspecs  o f  iflle  (first  tine  only) 


if (itime.eq.l)then 

string- ‘dir /si z/co:l/ou: transfer. tmp  • //if lie 
istat  us-1 lbs spawn (%descr (string (1 : len2 (string) )) ) 
if  ( .not . istat us l call  libSstop (Ival (1st at us) ) 


C  Read  garbage  from  the  directory  file 

open (unit -11, status- ‘old1 , file- ‘transfer. tmp' , fora-' formatted 
1  carriagecontrol-* list ‘ ) 

read  (11, •)  !  skip  empty  line 

read(H,*)  :  skip  dir  name  line 

read (11, •)  !  skip  empty  line 

end  if 

read (11, 10ClC,err-10) (string ( j : 3) , J-l, 21) 

10C1C  format  (2ial) 

if (string  (1 : 1) .eq. ‘  *)go  to  1C 
dev  flag  -  index  (if lie, ':* ) 
dc  3  -  1  ,  dev  flag 

wflle  ( 3 ;  J)  =  ifile(iO) 
end  do 

do  j  -  1, (index (string, ’;*) -1) 
k  -  5  ♦  dev_flag 
wflle(k:k)  -  strinQ(}:j) 
end  do 
ret  urn 

1C  stop  ‘No  more  files  match  the  input  spec.’ 

end 


B9  -  Routine  for  Data  Input  from  File 


Subroutine  NAME:  READER 


Written  by: 

Joseph  B.  Farrell 
DREA 

25  Aug.  1986 

Latest  revision:  3C  Aug.  1986 

subroutine  reader (blocka_to_read, bytes_per_biock,disk, raw_bytes, 
4  qult_f lag,  }rec, iblk, nblocks, quit_f iag2, lnput_channel,  ” 

4  disk~input_channel, blocka_per_aet ) 

IMPLICIT  INTEGER**  (a-*) 


byte 


raw_byte_data (50000) 


Integer*^ 
integer*4 
;nt eger#2 
integer *4 
integer*2 
integer*4 
integer*  4 
integer  *2 
integer *4 


blocks_tc_read 
blocks”t o”skip 
byte*_per_block 
disk 

input  channel 
disk  Tnput_channel 
raw_6yte§ 
text_losb(4) 
qult~f lag 
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integer**  quit_fiag2 

integer* 4  jrec 

integer**  iblk 

integer*4  nblocks 

common  /raw/raw_byte_data 

external  io$_skipf ile, lo$_readvblk, io$_skiprecord, ioS_rewind 


C —  Read  in  a  frame  of  data  from  tape  or  disk. 

if  ( jrec.eq.l)  parity_count  -  0 
if (disk. ne.l) then  “ 


Read  from  Magtape. 

do  i  -  1, block s_to_read 
parity_f lag  -  C~ 
mov  -  7i-l) ‘bytes  per_block  *  1 

status-sysSqiow (,Xval7input_channel) , io$_readvbik, 

1  text  iosb, , 

1  , %re7 (raw__byte_data (Mov) ) , %vai (bytes_per  Jdlock) , , , , ) 


C —  Check  to  see  if  there  was  an  error  on  the  tape  read. 

if (.not .status. or. text _iosb (2) ,eq.0)then 
frames_stored  -  jrec  -1 
call  set_cursor  (21, 1) 
if (text~iosb (4) ,ne.2.and.text_losb (45 
&  .ne. 10) then 

parity_count  -  parlty^count  ♦  1 


C—  If  parity  error  or.  the  first  record,  skip  into  the  file  until  we  get 
C —  tc  the  next  block  which  begins  with  the  first  channel,  then  begin. 

C 

if ( jrec.eq.l .and. i.le. block s_per_set) then 
type*, •  Error  reading  first  record:!* 
biocks_to_skip  -  biocks_per_set  -  i 
if (Diocks_tc_skip.ne.C) then 

status-sysSqiow  (, %val (input_channei) , 
i  io$  skiprecord,  text__iosb, , , 

4  %val (blocks_tc_skip7, , , , , ) 

if  (.not .status) call  lib$stop(%val (status) ) 
end  if 
go  to  92 
end  if 

type*,*  Error  reading  record::  Error  count-*, 

4  parity_count 


C—  If  a  parity  error  occurs  on  a  record  ,  go  back  and  reread  the  last 
C —  gooc  dIock  cf  data  which  begins  with  the  correct  channel. 


parity_flag  -  parlty_flag  ♦  1 
blocks_to_skip  -  - (fclocks_per  set  ♦  1) 
stat us-sysSqiow (, %val (input_cHannei) , 

4  io$  skiprecord, text_iosb, , , 

4  %vaT (biocks_to_skip) ,,,,,) 

if ( . not . stat us) call  1 ibSstop  (%val  (status) ) 
go  to  92 

else 

type*, 'End  of  file  encountered.* 
if  (raw_oytes . ne.C) then 
quit^flsQ  •  1 
return 
end  if 
end  if 

type*, 'Saving  * , f ramesstored, *  frames  and  exiting.* 
quit_fiag2  -  1 
ret  urn 
end  i  f 

raw_bytes  -  raw_bytes  ■*  bytes_per  clock 
i f  Tparity_f.ag.ne .Cither. 
bloc*s_tc  skip  •  parlty_fiag  •  b. ocks_per_set 
stat  us-sys'Jqlow  (, %val (input_channei ;  , 

4  io$  skiprecord, text_iosb, ,  7 

4  *vaT  (blocks_to_skipT, , , , , ) 

if  ( . net . status) cal  1  . IbSstop (%val (status) ) 
end  i  f 

end  ao 

(.St 


C-.etk  tc  see  if  the  next  read  will  put  us  past  the  end  cf  the  fi.e. 

is  tc  account  for  the  fact  that  DBREAE  counts  from  C  rather 
far  I.)  'Remember  that  we  haven’t  read  biock  iblk  yet.) 

if'' it. n*c.ocks_tc_read! . gt . (nblocks-1 ) ) then 


If  we're  out  cf  data,  dump  what  we've  collected,  or  slmp.y  close 
t  r.e  flies  if  we  have  no  data  tc  dunp, 

write  1  5,  *  *  Read  puts  us  pas*.  EOf.  b.ns  to  read  .  ‘ ,  pioe«cs  tc  read 
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write  (5,  • )  *  iblk  &  nblks  ‘ ,  lblk,  nblocks 
blocks_to  read  -  nblocks  -  iblk 
write  {5,  *) ’New^lks  to  read ‘ ,  Blocks_to_read 
pause 

type*, *  End  of  file  encountered.' 
if (blocks_to_read.gt .0) then 

raw_bytes  -  blocks_to_read  *  bytes_per_block 
quit_flag  “  1 
go  to  33 
else 

if(}.ne.l)  quit_flag2  -  2 
return 
end  if 
else 

C —  Read  from  Disk  file. 

C 

raw  bytes  -  blocks_to_read  *  bytes_per_block 
33  call  dbread (disk_input_channei, rawjoyte^dat a, blocks  tc_read, 

&  ierror,iblk)  - 

If  (ierror.ne.O) then 

typ>e*, 'Error  encountered  on  disk  read.' 
type* Saving  what  I  can  and  exiting.' 
qult_flag2  -  1 
return 
end  if 

lblk-iblk*blocks_to__read 
call  dbwait  <disk_input_channei) 
end  if 
end  i  f 


BIO  -  Transfer  Status  Routine 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

Subroutine  :  TRANSFER_STATUS 

Created  :  Summer,  1982  by  U.Vic  Physics  Co-op  student  Laurie  Bunch 

Major  Modifications  :  Spring,  1986  by  Joe  Farrell  -  Deep  Water  Acoustics 

Purpose  :  To  display  to  the  user  the  major  parameters  in  the 

Surveillance  Acoustics  TRANSFER  program  during  execution. 

Called  by  :  TRANSFER 

Calls  :  1)  DATE_TIM£  -  Subroutine  which  obtains  the  system  date 
and  time  in  an  ASCII  format. 

3)  ERASE_SCREEN  -  Subroutine  which  erases  the  screen  from 
the  specified  position  to  the  end. 

SET_CURSCR  -  Subroutine  which  sets  the  cursor  to  a 

specified  position  on  the  terminal  screen. 

5)  FORSSECNDS  -  Fortran  library  routine  which  aetermines 

the  number  of  seconds  difference  between 
the  number  specified  ana  tr.e  current  time. 

6)  COMPRESS  -  Subroutine  which  compresses  the  chosen 

channels  into  a  format  suitable  for  display. 


6)  COMPRESS 


Parameter  and  variable  definitions. 

REQUIRED  PARAMETERS: 

C KANN E  LS  _CH  0  S  EN  -  Channels  corresponding  to  the  hydrophones 
chosen  for  study. 

r.rec  •  Number  of  input  file  records  processed. 

FILES  •  Input  and  output  files. 

HYDROPHONES  CHOSEN  -  Hydrophones  chosen  for  study. 
NtJMBER_CHOSEn  -  Number  of  hydrophones  chcser  for*  study. 
numoer~cf_frames  -  Upper  limit  placed  by  tne  user  or.  the 
number  of  frames  to  be  procersec. 
f rames_procesed  -  Number  of  frames  which  have  teen  processea 
A LABEL-  •  ASCII  label  from  input  tape. 

7IME_OFFSET  •  Time  in  seconds  which  the  user  wishes  tc  step 
into  the  input  data. 

?C?AL_CHANNELS  »  Total  number  of  channels  or.  tr.«-  i -put  file. 
INTERNAL  VARIABLES: 

AVERAGE_TIME_per  frame  -  Amount  of  system  tir<  rot  CPU)  that 
one  frame  requires. 

BASE  -  The  number  of  rows  required  to  print  a.,  -a -or 
parameters  up  tc  the  channels  cnoser . 

BASE_7IME  •  Time  zero  when  the  timer  was  : - . * . a* r t . 
DATE-AND_7IM£  •  System  date  and  time  in  an  *S~  :  :  rmat . 

DELTA  TIME  -  Time  from  base  time  to  present. 

HMS_TTME_OFf SET  •  Hours,  minutes  and  second"  c?:  . .  v**er.t  of 
tne  time  offset  spec!'.*:  :  •  •  user  i ’ 

seconds . 
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FORMAT (/*  ' ,  T3,  ' Channel  usage  :  ‘,12,*  chosen  out  of  ' , 
12, 1  total1) 


Call  routine  to  compress  the  chosen  channels  into  a  format 
suitable  for  screen  output,  ie.  1-5, “7, 11, 15-24 . 


CALL  COMPRESS (Number  chosen, Channels  chosen, NELS) 
TYPE*,'  Channels  OUTVECT (:NEl3) 

CALL  COMPRESS (Number  chosen, Hydrophones_chosen, NELS) 
TYPE-,*  Acoustic  Channels:  \  OUTVECT { :NELS) 


Print  the  number  of  records  used,  the  number  of  FFT's  calculated, 
and  the  average  amount  of  system  time  per  FFT  per  channel  (not 
known  on  the  first  call.) 


TYPE  7000,  frames_processed 

FORMAT (/'  * ,T3, ‘Number  of  input  file  accesses  :  *,I6) 
TYPE  8000,  nrec 

FORMAT (/•  * ,T3, 'Number  of  records  processed  :  ',f6.2) 
TYPE  8100 

FORMAT ( '  ',T3, 'Average  time  per  frame  :  sec') 


Set  the  timer  and  find  the  absolute  value  of  time  zero. 


BASE  TIME  -  FORSSECNDS (0.0) 


For  successive  calls  just  print  the  dynamic  paramters. 


Find  the  amount  of  system,  time  to  have  elapsed  since  time  zero. 


DELTA  TIME  -  FOR5SECNDS (BASE  TIME) 


Calculate  the  average  time  per  FFT  per  channel. 


average_time_per_record  -  D£LTA_TIME  /  nrec 


Update  the  system  time. 


CALL  SET  CURSOR (2, 30) 
TYPE  9003,  DATE_AND  TIME 
FORMAT ('♦', A) 


If  the  number  of  channels  and  hydrophones  printed  on  the  first  call 
were  less  than  or  equal  to  16  then  16  lines  were  devoted  to  static 
parameters;  otherwise,  18  lines  were  devoted  to  static  parameters. 


IF  (NELS  .L£.  50)  THEN 
BASE  -  15 
ELSE 

BASE  -  17 
ENDIF 


Update  the  number  of  records  used,  the  FFT's  completed  and  the 
average  time  per  FFT  per  channel. 


CALL  SET  CURSOR (BASE  ♦  2,37) 

TYPE  9103,  frames_processed 
FORMAT ('♦', 16) 

CALL  SET  CURSOR (BASE  ♦  4,33) 

TYPE  9203,  nrec 
FORMAT ( ' ♦ ' , f  6 . 1 ) 

CALL  SET  CURSOR (BASE  +  5,27) 

TYPE  9303,  average  time  per  record 
FORMAT ( ' ♦ ' , F5. 2) 

ENDIF 


SUBROUTINE  COMPRESS (NUMBER_CHOSEN, INVECT, NELS) 


CHARACTER* 128  OUTVECT 


INTEGER*2 
INTEGER* 2 
INTEGER  2 
INTEGER *2 
INTEGER* 2 
INTEGER*2 
INTECER*2 
INTEGER*2 
INTEGER*2 


Number_chosen 

INVECT (NUMBER_CHOSEN) 

Last  flag 

START 

ST 

Temp  flag 
NELS 
L 
M 


COMMON  /OUT /OUTVECT 


1 f (number_chosen . eq . 1 ) t hen 

encodel2, 100,outveet (1:2) ) invect (1) 
format  (12) 

nels  -  2 
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r«t  urr. 

•  n<2  i  i 

START  -  INVECT  1 
ST  -  START 
NETS  -  C 

-  z 

Tamper. ag  -  Z 
K  -  2 

DC  J  •  2 ,  N-jmoar _cno»«r 

IF  (J.  ZZ.S'Jsxm:  cftonr.,^111  -  . 

IF 'last  r.ag.ES. : , THEN 

:f*:Rvect(J  .eg.ist*:.  then 

ST  -  ST  •  : 

CALL  BUILDER  1.  H.  NEDS,  I  .ag.  St  art .  5' 

ELSE 

TALL  BUILDER  'L.K.NELS.  T«*g  :„ag.  Start  St 
ST  -  INV ECTiJ 

CALL  BUI  LDER  NELS.  La«*  f.ag.St  .  St 

ENT  IF 
ELSE 

:r 'isr/Ecr  i:.  .es.  (ST*:  then 

ST  •  ST  •  : 

E15E 

TALI  BUILDER  :L.«,  NELS,  .^*t  :  .ag.  Start .  St 
START  -  IWVECT 
ST  -  START 
ENT  IF 
ENT  IF 
ENT  DC 
RETURN 
ENT 


SUBROUTINE  BUILDER  L,H. NELS, Laatf.ag, Start . S*  vp 

CHARACTER* 12*  OUTVECT 

CHARACTER*!  SEPARATOR 


INTEGER*: 
INTEGER *2 
INTEGER* 2 
INTEGER *2 
INTEGER*: 
INTEGER  *2 
INTEGER*: 


Last  * 
STAR* 
STOP 
Mcra  ? 
KELS’ 


•C 


*« 


CORMOK  /OCT /OUTVECT 

I F ( START . EC . STOP ) THEN 

ENCODE;:, ICC, OUTVECT ?L:Hi  STAR- 
FORMAT  (12; 

SEPARATOR  -  \  ’ 

CALL  FILLER (START, L.K. SEPARATOR,  LAST  FLAG.  NELS 
ELSE 

ENCODE ( 2 , ICC, OUTVECT (L:M  STAR- 
SEPARATOR  - 

CALL  FILLER (START, L.K, SEPARATOR.  LAST  FLAC.NEIS 
encode cutvect .l:k  . stop 

SEPARATOR  - 

CALL  FILLER 'STOP, L,K, SEPARATOR.  LAI"  r  LAC.  HE  _S 
END  IF 
RETURN' 

END 


SUBROUTINE  FI  LLER  1  ELEMENT  ,  L,  *.  SEP  AAATCR  .  -AS’  t  1>G 


CHARACTER  *128  OUTVECT 

CHARACTER*:  SEPARATOR 


INTEGERS 
INTEGER*: 
INTEGER *2 
INTEGER*: 
INTEGER*: 


La at  t. 

ELEMENT 

NELS 


•S 


k 


C  OWON  t  OUT  /  OUTVECT 

IF  f ELEMEKT.LT.::; THEN 

CUTVECT  ' L :  Ll  -  CUrVEC” fM:K 
SFfLatt  .AND. SEPARATOR. EL  "Mf* 

NELS'-  NELS  *  : 

ELSE 

OUTVECT (K:K  -  SEPARATOR 
NELS  -  NELS  *  2 

k  •  k  ♦  : 

END  IF 
ELSE 

IF  (La»*._f.ag.  EC.  I  .AND  .  SEPARATOR  i; 

NELS  -  NELS  •  2 
ELSE 

OUTVECT  •  (M*: »:  (M-: .  -  SLPAfiATDR 


.  ne: 
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